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FIXED-BASE SIMULATION EVALUATION OF VARIOUS 
LOW-VISIBILITY LANDING SYSTEMS FOR HELICOPTERS 

By P.S. Rempfer, L.E. Stevenson, and J.S. Koziol, Jr. 
Electronics Research Center 


SUMMARY 

A fixed-base simulation evaluation of various low-visibility 
landing systems for helicopters has been made. The low-visibility 
mission consisted of a straight-in glide slope and localizer 
approach with a flare and deceleration maneuver at the end. The 
landing systems consisted of a fully automatic system and six 
manual systems. The manual systems consisted of three flight 
control modes each being flown with a flight director indicator 
and then with raw data displays. The three flight control modes 
were an attitude command mode with an unaugmented vertical axis, 
an attitude command mode with rate of descent command in the 
vertical axis, and a three-axis velocity command mode. The 
landing systems and the helicopter were simulated on a hybrid 
computer. The landing systems were al] digital. A series of six 
subjects consisting of two engineers and four instrument-rated 
pilots flew each of the modes for the landing mission. Per- 
formance data was recorded by the digital computer and analog 
recorders. The relative performance of the modes are presented. 

The results indicate that all modes are acceptable for the 
mission prior to flare and deceleration. When flare and decel- 
eration are considered, the simplest attitude command mode 
becomes unacceptable both with a flight director indicator and 
with raw data displayed. This is due primarily to the fact that 
this mode had no augmentation in the vertical axis and the 
coupled dynamics between pitch and power became difficult to 
control through the flare and deceleration maneuver. When 
augmentation in the vertical axis was added to provide a rate of 
descent command system with altitude hold capability, the system 
became acceptable both with the flight director indicator and 
with raw data displays. 


INTRODUCTION 

The majority of helicopters flying today do so with no 
electronic controls or, at most, with limited authority attitude 
rate stabilizing systems. Although these vehicles demand a high 
pilot workload they have been very successfully used under VFR 
conditions and under IFR conditions in cruise flight. Low-speed 
flight under IFR has been carried out by the military to some 
degree but it is generally agreed that advanced systems in con- 


trols, displays, and navigation aids are required before the 
full range of IFR helicopter operation can be realized. 

Recognizing this, research and development of advanced IFR 
helicopter systems was initiated about ten years ago by the 
military and civilian agencies in the U.S. and abroad. The 
approach taken by the different agencies varied, resulting in a 
variety of control/display systems. For example, Langley 
Research Center has done work on a variety of displays with 
attitude rate stabilized vehicles (ref. 1) and on full authority 
attitude command control with flight director indicator displays 
(ref. 2) . The British Royal Aircraft Establishment has studied 
flight director display with attitude rate stabilized vehicles 
(ref. 3). The Air Force Flight Dynamics Laboratory has worked 
with an attitude command control system with augmentation in 
the vertical axis to give rate of descent command/attitude hold 
(ref. 4). The U.S. Army Transportation Research Command is also 
considering a three-axis translational velocity command control 
system (ref. 5). Before commercial helicopter operations are 
able to include low-speed low-visibility missions, one of these 
types of advocated systems must be selected. The criteria 
for selection will involve trade-off studies among cost, system 
reliability, mission reliability, and most certainly system 
performance . 

The objective of this investigation is to present performance 
data on a number of low-visibility landing systems to be used 
in making these required trade-off studies. The low-visibility 
landing mission which was considered is shown in Figure 1. 

A conventional localizer intercept is made at approximately 800 
ft of altitude and the helicopter is decelerated to 42 knots 
range rate. The glide path is intercepted and followed at 42 
knots range rate. At 150 ft of altitude a combined flare and 
deceleration is made to 50 ft of altitude and 10 knots range rate. 
This condition is held while tracking the localizer until at 100 
ft of range-to-go visual contact is assumed to have occurred. 

In order to study a number of landing system concepts, a 
flexible digital multi-mode landing system was developed. This 
system allowed the pilot to select any one of a variety of 
landing systems. These landing systems consisted of a fully auto- 
matic system and six manual systems. The manual systems consisted 
of three flight control modes each being flown with a flight 
director indicator and then with raw data displays. The three 
control modes were an attitude command mode with unaugmented 
vertical axis, an attitude command mode with rate of descent 
command and a three-axis translational velocity command mode. 

The helicopter along with the landing system was simulated 
on a hybrid computer. The helicopter simulated was the CH-46C 
helicopter stationed at NASA Langley Research Center. This 
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vehicle was chosen since it was available for possible flight 
test of the landing system. The pilot interface was provided with 
a single seat fixed base cockpit. The digital multi-mode landing 
system was simulated on the digital portion of the hybrid com- 
puter. The functional layout of the simulation is shown in 
Figure 2. 



The fixed-base simulator studies were conducted using both 
instrument-rated pilots and engineers with simulator experience. 
The low-visibility mission was flown several times for each 
landing system by each subject. The resulting data are presented 
in this report along with some preliminary conclusions. 
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Figure 2.- Hybrid computer simulation 


The report is organized such that the reader may get an over- 
view by reading the Introduction, Results, and Conclusions. The 
details of the low-visibility mission and its selection are pre- 
sented in the Low-Visibility Mission Section. The details of 
the multi-mode landing system and a tutorial discussion of its 
synthesis are presented in the System Description Section. The 
details of the hybrid simulation are presented in the Simulation 
Description Section. Finally, the experimental procedures used 
in making the tests and obtaining the data are detailed in the 
Experiment Section. 

These studies were conducted at the Electronics Research 
Center in the period March, 1968 to January, 1970. Although the 
authors were the principal investigators , many people assisted in 
carrying out the effort. In addition to the NASA employees, 
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involved in these studies were the following: Messrs. Stanley 

Driban, Mannie Smith, Jay Haas, and Roger Brown of Service 
Technology Corporation were responsible for the hybrid programming 
and cockpit interface throughout the effort; mssrs. Joseph P. 
Tymczyszyn, Thomas Imrich, graduate students at M.I.T., and 
John Spencer of Bell Aerosystems acted as test subjects. The 
authors gratefully acknowledge the assistance of these people 
without whose help the effort would not have been possible. 


SYMBOLS 

A , A , A , A Actuator output state in the vertical, 

pitch, roll, and yaw axes, respectively, 
inches 

DCC, DEC, DAC , DRC Control system input in the vertical, 

pitch, roll, and yaw axes, respectively, 
inches 


DCP , DEP, DAP, DRP Pilot inputs from collective stick, 

longitudinal center stick motion, lateral 
center stick motion, and pedals, respectively, 
inches 


DCCBIAS 

DCR 

DCRO 

DER 

DERO 

DV2 

DV4 

DV5 

DV10 

E 

ENOM 

EMAX 


Reference control input in vertical axis, inches 
Reference collective stick input, inches 
Reference collective stick input, inches 
Reference longitudinal stick input, inches 
Reference longitudinal stick input, inches 
Velocity command limit, ft/sec 

Forward loop gain in altitude loop of guidance, sec“l 
Velocity command rate limit, ft/sec^ 

Forward loop gain in altitude loop of FDI , sec ^ 
Variable in performance index 

Nominal value of a variable in the performance index 

Maximum allowable excursion of a variable in the 
performance index from its nominal value 


5 



g 

Gg 

G$ 

GW 

GQ 

GR 

GS 

G0 

GY 

GZ 

GI$ 

GIG 

GVXI 

GVX 

GVY 

GZD 

g xa' 

GV1 


Acceleration of gravity at surface of the earth, 
ft/s ec^ 

Forward loop gain in beta loop of yaw axis, sec ^ 

Forward loop gain in roll loop of roll axis, 
in/radian 

Forward loop gain in heading loop of yaw axis, sec ^ 
Feedback gain on pitch rate, sec 

Forward loop gain in yaw rate loop of yaw axis, in/ 
(radian/sec) 

Glide slope angle, radians 

Forward loop gain in pitch loop of pitch axis, 
in/radian 

Forward loop gain in lateral position loop of roll 
axis, sec _ l 

Forward loop gain for altitude loop of vettical 
axis, sec - ! 

Integral-by-pass gain in roll loop of roll axis, 
sec-1 

Integral-by-pass gain in pitch loop of pitch axis, 
sec - ! 

Integral-by-pass gain in velocity loop of pitch axis, 
radians/ft 

Forward loop gain in velocity loop of pitch axis, 
radians/ (f t/sec) 

Forward loop gain in velocity loop of roll axis, 
radians/ (ft/sec) 

Forward loop gain in velocity loop of vertical axis, 
in/ (ft/ sec) 

GyA , G za Wind gusts along the X, Y, Z, body axes, 
respectively, ft/sec 

Inteqral-by-pass gain in altitude loop of guidance, 
sec - l 
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GV2 


GCFR 

GCF'!' 

GLAATT 

GLOATT 

GVATT 

GWEL 

1 , 1 , 
xx ' yy ' 

I 

xy 

J 

xy 

L, M, N 
m 

P, Q, R 

RFDI 

RFDT 

S 

S$ 

Sty 

SZD 

SVX 


Forward loop gain in altitude loop of guidance, sec 

Roll command to yaw crossfeed in yaw axis, sec 

Roll command to yaw crossfeed in yaw axis of automatic 
mode, sec"l 

FDI input gain in roll loop of roll axis, volts/radian 

FDI input gain in pitch loop of pitch axis, 
volts/radian 

FDI input gain in raw collective loop of vertical 
axis, volts/in 

FDI input gain in velocity loop of vertical axis, 
volts/ (ft/sec) 

I Moment of inertial about the X, Y, Z, body axes, 
respectively, slug ft^ 

Product of inertia = J ' xy dm, slug ft^ 

Product of inertia = -J" xy dm, slug ft^ 

Applied moments about the X, Y, Z, body axes, 
respectively, ft-lbs 

Mass of the helicopter, slugs 

Angular rate of helicopter along the X, Y, Z, body 
axes, respectively, radians/sec 

Range at initiation of deceleration, ft 

Range at termination of deceleration, ft 

Laplace transform variable 

Stick sensitivity in roll loop of roll axis, radians/in 

Pedal sensitivity in heading loop of yaw axis, 
(radians/sec) /in 

Collective stick sensitivity in velocity loop of 
vertical axis, (ft/sec) /in 

Stick sensitivity in velocity loop of pitch axis, 
(ft/sec) /in 
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SVY 


S0 

SB 

TAS 

TAST 

TB 

t 

TAU 

TVZD 

TCFR 

TCFY 

TP 

U, V, W 
VXCR 

VXCRO 

VXCI 

VXCF 


Stick sensitivity in course loop of roll axis, 
(radian/sec) / in 

Stick sensitivity in pitch loop of pitch axis, 
radian/in 

Pedal sensitivity in beta loop of yaw axis, radian/in 
Airspeed of the helicopter, ft/sec 

Airspeed at which transition from turn coordination 
to heading hold in yaw axis takes place, ft/sec 

Lag time constant in beta loop of yaw axis , sec 

Current time in wind noise generation, sec 

Elapsed time since flare initiation, sec 

Total time of open loop flare, sec 

Lag time constant on roll command to yaw crossfeed, 
sec 

Lag time constant in roll command to yaw crossfeed in 
automatic system, sec 

Roll rate feedback gain, sec 

Inertial velocity of the helicopter along the X, Y, Z, 
ft/sec 

Reference forward velocity in velocity loop of the 
pitch axis, ft/sec 

Reference forward velocity in velocity loop of FDI , 
ft/sec 

Velocity command at initiation of deceleration, ft/sec 
Velocity command at termination of deceleration , ft/sec 


V , V , V Inertial velocity of helicopter along the X, Y, Z, 
x Y z axes of the approach navigation frame, respectively, 
ft/sec 



V h Inertial velocity of helicopter along the X, Y, Z, 
z axes of the heading vertical frame, respectively, 
ft/sec 


VLAI FDI input to vertical command bar, volts 

VEI FDI input to vertical tab, volts 
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VLOI 


FDI input to horizontal command bar, volts 


VXCO Initial velocity command in vertical axis for flare, 

ft/sec 

VZGR Vertical velocity reference in guidance, ft/sec 

X, Y, Z Position of the helicopter in the approach navigation 
frame when not subscripted, feet. When subscripted, 
these variables represent the applied forces along 
the X, Y, Z body axes, respectively, lbs 

XFLI Range at initiation of flare, feet 

ZFLT Altitude at initiation of flare, feet 

ZR Reference altitude in altitude loop of vertical 

axis, feet 

ZRO Reference altitude in altitude loop of FDI, feet 

ZCGS Glide slope altitude at current range, feet 

a n Correlation time of noise process 

3 Sideslip angle, radians 

3 n Noise process parameter 

Y Locke number 

T Noise process variable 

6 , 6 , 6 , 6 Rotor output variables in the vertical, pitch, 

roll and yaw axes, respectively, inches 

A n Noise sequence discretization interval 

C Actuator model damping ratio 

Cl 

C n White noise with spectral density of unity 

n n White quassian random sequence with zero mean and 

unity variance 

¥, 0, $ Yaw, pitch, and roll Euler angles taken in that order 
from a local vertical frame to the helicopter body 
axes, radians 

0R Reference pitch attitude, radians 
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G 

X 

CO 

a 

A ( ) 
VAR ( ) 
/ ( ) 

Z ( ) 

I ( ) I 

E ( ) 

d_( ) 
dt( ) 

i (*) 

9 ( ) 

m 

( ) 


Course of helicopter, radians 

Reference course of helicopter in course loop of roll 
axis, radians 

Noise process variable 

Variance of noise process 

Noise process parameter 

Dummy variable of integration 

Correlation time 

Reference heading for heading loop of yaw axis, radians 
Actuator model natural frequency, rad/sec 
Angular rate of the rotor, rad/sec 
Perturbation of ( ) 

Variance of ( ) 

Integral of ( ) 

Summation of ( ) 

Absolute value of ( ) 

Expected value of ( ) 

Total time derivative of ( ) 

Partial derivative of (*) with respect to ( ) 

Total time derivative of ( ) 


Subscripts 

p, q, r, u, v, w, 6 subscripts on X, Y, Z,; L, M, N, denote the 

partial derivative of the aerodynamic force 
of moment with respect to P, Q, R, U, V, W, 
<$ c , S e , S a , or <S r , respectively, (i.e., 
stability and control derivatives) . 
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A subscript on X, Y, Z, denotes the force is due only 

to aerodynamics 

TR denotes the variable is a trim value in steady level 

flight 

COM denotes the variable is a command input 

IC denotes initial value or condition 

AS denotes velocity with respect to air mass rather 

than inertial 


Axis Systems 

Body Frame . - Origin is fixed as the normal center of gravity 
of the helicopter; X-axis is directed forward out the nose of the 
helicopter, and parallel to the water line; Z-axis is directed 
downward and perpendicular to the water line; Y-axis completes 
the orthogonal triad. 

Approac h Nav igation Frame .- Origin is fixed at the runway 
touchdown point; X-axis is horizontal and directed in the direc- 
tion of travel for an approach; Z-axis is vertical and directed 
downward; Y-axis completes the orthogonal triad. 

Vert ical Heading F r ame .- Origin is fixed at the normal center 
of gravity of the helicopter; X-axis is the orthogonal projection 
of the body X-axis onto the local horizontal plane; Z-axis is 
vertical and directed downward; Y-axis completes the orthogonal 
triad . 


LOW-VISIBILITY MISSION 

The low-visibility mission which was considered consists of 
a straight-in approach and landing. The mission profile appears 
in Figure 1. The mission assumes visual contact with the desired 
touchdown point will occur before 100 ft of range-to-go at 50 ft 
of altitude. With landing lights available, this is virtually a 
zero visibility condition. A conventional localizer intercept is 
made at 800-ft of altitude. A 0.1 radian (6-degree) glide path 
is intercepted and tracked at 42-knots range rate. At 150 ft of 
altitude, a flare and deceleration maneuver is made to 50 ft 
altitude and 10 knots range rate. This condition is held with 
localizer track until visual contact is made. 

This mission is conservative in that it is a simple extension 
of current fixed wing aircraft IFR operation. It requires 
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approximately 1500 ft of range without obstructions protruding 
above the height of the desired touchdown point. Such a mission 
could be flown into ports over a body of water, to the roofs of 
high buildings or, possibly, to ports built over freeways with 
the approach made along the freeway. Depending on the exact 
nature of the mission, either ILS-type equipment or approach radar 
could be used. Since the ILS-type equipment is the most con- 
straining, it was chosen as the basic approach navigation source 
for this study. The profiles generated are compatible with a 
radar source of information. 

The parameters of the mission trajectory were selected with 
the following considerations: 

(a) A glide path of 0.1 radian was used based upon the 
material presented in reference 6. The glide path 
source is located 50 ft in front of the desired touch- 
down point so that the aircraft will fly over the 
source making it available for the next aircraft. 

During the simulation, the conventional fixed wing 
glide path sensitivity of ± .5 degrees full-scale was 
found to be too sensitive for this mission. A 
sensitivity of ± 2.0 degrees full-scale was found by 
experiment to be acceptable to within 1500 ft of range. 

(b) The range rate along the glide path is 42 knots so that 
the aircraft will not exceed the 500 to 700 fpm rate of 
descent limit sited in reference 7. This range rate is 
also in keeping with the more recent studies of 
reference 2. 

(c) The flare and deceleration are initiated at 150 ft of 
altitude. This altitude was selected based upon a 
tradeoff between the following: 

(1) Hold maximum range rate as long as possible to 
get down quickly, 

(2) Keep the range over which a radar altimeter might 
have to be used to a minimum, 

(3) Stay within the 0.15-g deceleration limit recommended 
in reference 8. 

(4) Keep any dipping below the desired final altitude 
to a minimum. 

(d) The final altitude is 50 ft which is consistent with 
reference 2. This altitude might have to be adjusted 
for operational use. 
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(e) The final range rate was chosen at 10 knots. This is 
a modest rate from which visual takeover at 100 ft of 
range could be smoothly accomplished. 

(f) The glide slope intercept altitude is 800 ft. This 
corresponds to a range of 8050 ft. From intercept to 
flare initiation is 6500 ft. At 42 knots, this gives 
the pilot just over the 1.5 minutes recommended in 
reference 6 to stabilize on the glide path before 
flare . 

(g) The localizer source is located 1000 ft down range. 
During the simulation, the conventional fixed wing 
localizer sensitivity of ± 2.5 degrees full-scale was 
also found to be too sensitive for this mission. A 
sensitivity of ± 5.0 degrees full-scale was found by 
experiment to be acceptable. The combination of the 
resultant localizer sensitivity with the location of 
the localizer source gave full scale over the desired 
touchdown point of ± 100 ft assuring a pilot on scale 
of visual contact. At 10,000 ft-range, the localizer 
behaves very much like a standard ± 2.5-degroe 
sensitivity localizer at the far end of a 10,000-ft 
runway . 

(h) The range and range rate source are located with 
the localizer for convenience. 


SYSTEM DESCRIPTION 

In this section, the digital multi-mode landing system modes 
are defined; the block diagrams and gain values for each mode are 
presented, and the flight control system synthesis based upon 
linear analysis and simulation is shown. The flight director 
laws were synthesized with reference to the automatic mode. 


Mode Definition 

The digital flight control system was comprised of four 
flight control (FC) modes and three associated flight director 
indicator (FDI ) modes. Each FC mode was a four-axis system 
controlling differential collective rotor pitch (DEC) for pitching 
moment, collective rotor pitch (DCC) for lift, lateral cyclic 
rotor pitch (DAC) for rolling moment, and lateral differential 
cyclic rotor pitch (DRC) for yawing moment. The DRC input is 
analagous to rudder input in fixed wing aircraft. Electro- 
hydraulic actuators of an Electric Input Servo System (EISS) 
served to force the required control linkages. Each FC mode 


13 



accepted four axes of pilot input. Each input device produced 
an electric signal and had artificial force feel. The inputs 
were longitudinal motion of center stick (DEP) normally to DEC, 
lateral motion of center stick (DAP) normally to DAC , pedal 
motion (DRP) normally to DRC, and collective stick motion (DCP) 
normally to DCC . The yaw axis for the FC modes varied with air- 
speed having a low-speed characteristic and a high-speed charac- 
teristic . 

The FDI is shown in Figure 3. Each FDI mode was a three- 
axis system addressing the horizontal command bar (VLOI) for 
pitching motion commands, the vertical "tab" (VEI) for lift 
commands and the vertical command bar (VLAI) for rolling motion 
commands . 



Figure 3.- Flight director indicator 
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The sense of the FDI commands, the pilot inputs, and EISS 
inputs are defined in the following tables. 

PILOT INPUTS 


DCP 

gradient 

is 

positive 

for 

stick 

down 

DEP 

gradient 

is 

positive 

for 

stick 

back 

DRP 

gradient 

is 

positive 

for 

right 

pedal in 

DAP 

gradient 

is 

positive 

for 

stick 

right 


EISS INPUTS 

DCC positive gives an upward acceleration 
DEC positive gives a pitch up acceleration 
DRC positive gives a yaw right acceleration 
DAC positive gives a roll right acceleration 

FDI INPUTS 

VEI positive moves the vertical tab down 
VLOI positive moves horizontal bar down 
VLAI positive moves vertical bar right 


A definition of each of the FC modes and associated FDI 
modes follows, along with their block diagrams. These definitions 
were arrived at by considering the complexity of the sensors and 
computations which would be required to implement each mode in 
operation. Hypothesized sensors for each mode are presented in 
Appendix A. 

In each manual mode, the pilots' command inputs from the 
center sticks and the pedals are passed through a dead zone to 
eliminate small amplitude noise and bias inputs. The dead-zone 
transfer function is shown in Figure 4. It is essentially a 
unity pass transfer function with threshold. It should be noted 
that added variables have been introduced into the various modes 
and given values when the particular mode is first entered. These 
references are implemented to prevent system transients when the 
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desired mode is initially entered or changed. Finally, in each 
of the flight control modes and flight director modes all inte- 
grators are initially set to zero when the particular mode is 
first entered. 



Figure 4.- Dead zone transfer function 


Attitude Command Mode I.- This mode is basically an attitude 
command mode Tn roll and in pitch. The vertical axis is unaugmented 
and the rudder axis is a sideslip command system at high airspeeds 
and a heading rate command system at low airspeeds. The block 
diagrams are given in Figures 5—11. A precise definition of the 
control and FDI inputs follows. 

DEP input is interpreted as a pitch attitude command. 

The helicopter responds by adopting and holding a 
pitch attitude proportional to the DEP input. 

DAP input is interpreted as a roll attitude command. 

The helicopter responds by adopting and holding a 
roll attitude proportional to the DAP input. 

DCP input is interpreted as a raw collective rotor pitch 

input. The helicopter is not augmented in this axis. 

DRP input interpretation varies with airspeed. At low 

airspeeds the input is interpreted as a heading rate 
command. The helicopter responds by tracking a 
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GQ = 0.5 SEC. 

G® = 20. IN/RADIAN 

®R S (®+DEC/G®)(AT MODE ENGAGE) GI® = 0.2 SEC.” 1 

DER = DEP (AT MODE ENGAGE) S® = 0.1 RADIAN/IN 


Figure 6.- Pitch axis of Attitude Command Mode I 
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T/3 = 0.5 SEC 

S/3 = -0.2 RAD/IN 

6/3 = 0.3 SEC -1 

S<X> = 0.1 RAD/IN 

GCFR = 0.3 SEC" 1 
TCFR = 0.5 SEC 
GR = 15. IN/(RAD/SEC) 

Figure 7a.- Yaw axis of Attitude Command Mode I (TAS > TAST) 



¥ R =¥(AT MODE ENGAGE) 

S¥ = 0.2 (RAD/SEC)/ IN 

GV = 1.0 SEC' 1 

GR = 15. IN/(RAD/SEC) 

TAST = 50. FT/SEC 

Figure 7b.- Yaw axis of Attitude Command Mode I (TAS < TAST) 
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S<D = 0.1 RAD/IN 
TP = 0.5 SEC 
Gtt> = 0.35 SEC" 1 
G<X> = 15. IN/RAD 


Figure 8.- Roll axis of Attitude Command Mode I 



DERO * DEP (AT MODE ENGAGE) 
PITCH RATE LIMIT * ±0.1 RAD/SEC 
VLOI FULL SCALE =±2.6 VOLTS 


Figure 9.- VLOI/FDI axis of Attitude Command Mode I 
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VLAI FULL SCALE* ±2.6 VOLTS 



VLAI 


IDAP 


GY *0.2 SEC' 1 
GVY * .04 RAD /(FT/SEC) 
S<S> = .1 RAD/IN 
GLAATT* 3. VOLT/RAD 


Figure 10.- VLAI/FDI axis of Attitude Command Mode I 



ZRO = Z (AT MODE ENGAGE) 
DCRO = DCP (AT MODE ENGAGE) 

VEI FULL SCALE = ±16. VOLTS 


DVIO = 0.3 SEC' 1 
GZD = -0.2 I N/( FT/SEC) 
GVATT = -5. VOLT/ IN 


Figure 11.- VEI/FDI axis of Attitude Command Mode I 
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heading rate proportional to the DRP input. Zero 
input is heading hold. At high airspeed, the input 
is interpreted as a sideslip command. The helicopter 
responds by adopting and holding a sideslip angle 
proportional to the DRP input. Zero input provides 
for a coordinated turn. The airspeed value separating 
low and high speed is called the transition airspeed 
(TAST) . 

VLOI output represents the error between the pitch 

attitude desired by the landing system outer loops 
and the current pilot-commanded pitch attitude. 

VLAI output represents the error between the roll attitude 
desired by the landing system outer loops and the 
current pilot-commanded roll attitude. 

VEI output represents the error between the collective 

rotor pitch desired by the landing system outer 
loops and the current pilot-commanded collective 
rotor pitch. 


Attitude Command Mode II .- This mode is similar to the 
Attitude Command Mode I except that the vertical axis is a rate 
of descent command system with altitude hold capability. The 
block diagrams for the vertical axis are given in Figures 12 
and 13. A precise definition of the control and FDI inputs for 
the vertical axis follows. 

DCP input is interpreted as a rate of descent command. 

The helicopter adopts and holds a rate of descent 
proportional to the DCP input. Zero input provides 
an altitude hold. 

VEI output represents the error between the rate of 

descent desired by the landing system outer loops 
and the current pilot-commanded rate of descent. 

Velocity Command Mode .- In this mode the pilot controls his 
translational velocity with respect to the ground. The vertical 
axis and rudder axis are the same as in Attitude Command Mode II. 
Forward displacement of the center stick results in rate of 
change of the helicopters forward velocity. Lateral displacement 
of the center stick results in a rate of change of the helicopter's 
course along the ground. Due to the nature of the rudder axis 
this course change is accomplished with a coordinated turn at 
high airspeed and a sideslip maneuver at low airspeed. The block 
diagrams required by this mode are presented in Figures 14 to 17. 
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SZD=8.(FT/SEC)/IN 

GZ = 1.3 SEC -1 

GZD = -.2 I N/( FT/SEC) 

DCCBIAS = DCC (AT MODE ENGAGE) 

ZR = Z (AT MODE ENGAGE) 

DCR= (DCP-V z /SZD)(AT MODE ENGAGE) 

OUTPUT OF THE ALTITUDE INTEGRATOR RESET 
TO -150. FT AT FLARE INITIATION 

Figure 12.- Vertical axis of Attitude Command Mode II 



DCRO = (DCP-VZ/SZD)(AT MODE ENGAGE) 
VEI FULL SCALE = ±I6.V0LTS 


Figure 13.- VEI/FDI axis of Attitude Command Mode II 
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VXCRMVx COS'P+Vy SIN 't )( AT MODE ENGAGE) 
DER * DEP (AT MODE ENGAGE) 

OR* (8 + DEC/G8MAT MODE ENGAGE) 

RESOLVER •• 

Vx = V x COS¥ + V Y SIN ¥ 


SVX--3.0 (FT/SEC 2 )/IN 

GVX* -.01 RAD/ (FT/ SEC) 

GVXI* -.001 SEC -1 

GQ * 0.5 SEC 

G8- 20. IN/ RAD 

GI8 = 0.2 SEC' 1 

S8* 0.1 RAD/ IN 


Figure 14.- Pitch axis of Velocity Command Mode 



SVY < 


.04 RA °' S -~ i v *" t TO. FT / SEC 

.0133 <2.+ 70./vJ ) RA P /SEC ; IO.FT/SEC«vJ <70.FT/SEC 
* IN 

0 RAO/SEC . v h s ,o. FT /SEC 
IN 

RESOLVER' 

V* ■ (V, COS ¥+ V y sin¥) 

v? ■ ( v„ cos'i'-v. siny ) 

y j " 


S<1> = 0.1 RAO/ IN 
GVY s .04 RAO/ (FT/ SEC) 
TP*O.S SEC 
GI$*<X35 SEC' 1 
G<X><I5.0 IN/RAD 


Figure 15.- Roll axis of Velocity Command Mode 
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DERO > DEP (AT MODE ENGAGE) 

VXCRO-(V x COS¥+V y SIN''IO<AT MODE ENGAGE) 

VLOI FULL SCALE * ± 2.6 VOLTS 
RESOLVER 

vj- V„COS¥+ V y SIN^ 
vj- V y COS^- V X SIN^ 

Figure 16.- VLOI/FDI axis of Velocity Command Mode 


GVX * - .01 RAD/IFT/SEC) 
GVXI *-.001 RAD/FT 
SVX— 3.0(FT/SEC Z )/IN 
Se* 0.1 RAD/IN 
GLOATT- -5. VOLT/RAD 
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(^t-gy -k^t- gvy — [ 


6LAATT 


VLAI 


DAP 


DEAD 

ZONE 


■EH 


SVY - 7 V** TAN (£-¥> -H^-JCVY 

V tt— 

ifT~L 


GVY [-Hgh 


£R * (TAN ' ( V y h /V x h ) + ¥)(AT MODE ENGAGE) 


RESOLVER 


SVY = 


04 M2£SEC . v h i 70>FT/SEC 


S4> *0.1 RAD/ IN 


.0133 <2.+ 70./V„ h } RA °^, S — ; IO.FT/SEC<V„ h < 70. FT/SEC GVY *.04 RAD /(FT/SEC) 

X ,N GY* 02 SEC' 1 

0 RAD/SEC . yh s | 0 FX/SEC GLAATT * 3.0 VOLT/RAD 

I N * 

VLAI FULL SCALE * - 2.6 VOLTS 


RESOLVER 

v x * ( VjCOS'I' + VySin'I') 

v y h * ( v y cos 'J'-v* sin'll 


Figure 17.- VLAI/FDI axis of Velocity Command Mode 

A precise definition of the control and FDI inputs for this 
mode follows. 

DEP input is interpreted as a rate of change command 

on the longitudinal translational velocity of the 
aircraft along the X-axis of the vertical heading 
axis system. The helicopter responds to the input 
by adopting and holding a rate of change of this 
velocity proportional to the DEP input. Zero 
input is a velocity hold. 

DAP input is interpreted as an angular rate of change 

command on the course of the aircraft. The helicopter 
responds to the input by adopting and holding a rate 
of change of course proportional to the DAP input. 

Zero input is a course hold. 

VLOI output represents the error between the pitch attitude 
desired by the landing system outer loops and the 
pitch attitude currently being commanded by the 
longitudinal velocity rate command system. 
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VLAI 


output represents the error between the roll attitude 
desired by the landing system outer loops and the 
roll attitude currently being commanded by the course 
rate command system. 

Automatic Mode.- The automatic mode does not accept pilot 
command signals. The FDI is inoperative with the pilot monitoring 
the mission with raw data displayed via conventional instruments . 
The instrument panel layout is shown in Figure 18. The block 
diagrams are presented in Figures 19 to 22. The longitudinal and 
vertical guidance laws required by this mode and the FDI modes are 
presented in Figures 23 and 24. 





'GLIDE SLOPE’ 
AND 
LOCALIZER 
DEVIATION 
INDICATOR 


DIGITAL 

RANGE 

INDICATOR 


Figure 18.- Instrument panel layout 
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ZR = Z (AT MODE ENGAGE) 

DCCBIAS = ( DCC + (GZD) V Z )(AT MODE ENGAGE) 


Figure 19.- Vertical axis of Automatic Mode 



RESOLVER « 

vj = V x COS ¥ + V y SIN ¥ 

PITCH RATE LIMIT 5 ±0.1 RAD/SEC 
0R = (8+ DEC/G8MAT MODE ENGAGE) 


GVX = -.01 RAD/(FT/SEC) 
GVXI = -.001 RAD/FT 
GQ = 0.5 SEC 
G0 = 20. IN/RAD 
GI0 = 0.2 SEC -1 


Figure 20.- Pitch axis of Automatic Mode 
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ROLL RATE LIMIT = ±0.5 RAD/SEC 


GY = 0.2 SEC -1 

GVY = 0.04 RAD/( FT/SEC) 

TP = 0.5 SEC 

GI<X> = 0.35 SEC -1 

G4> = 15. IN/RAD 


Figure 21 .- Roll axis of Automatic Mode 
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GVY = 0.04 RAD/(FT/SEC) 

ROLL RATE LIMIT = ±0.5 RAD/SEC GCFR = 0.3 SEC -1 


TCFR = 0.5 SEC 


GR = 15. IN/(RAD/SEC) 


TAST = 50. FT/ SEC 


Figure 22a.- Yaw axis of Automatic Mode (TAS > TAS.T ) 



ROLL RATE LIMIT ■ ±0.5 RAD/SEC G ' i ' * 10 < RAD/SEC > /RA0 

GR* 15. IN/(RAD/SEC) 


YR* Y (AT MODE ENGAGE) 


TAST * 50. FT/SEC 


Figure 22b.- Yaw axis of Automatic Mode (TAS < TAST) 
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[vxci-vxcf][i-3 (r P^fdt ) 2 + 2 (iSr) 3 ] + Vf 


VXCI 

LONGITUDINAL 
VELOCITY 
COMMAND 

VXCF 

RFDI RFDT 0 

RANGE 
(X) 

RFDI = -1550. FT. (RANGE AT INITIATION OF DECELERATION) 

VXCI = 70.FT/SEC (VELOCITY COMMAND AT INITIATION OF DECELERATION) 
RFDT? -275. FT. (RANGE AT TERMINATION OF DECELERATION) 

VXCF = 15. FT/SEC (VELOCITY COMMAND AT TERMINATION OF DECELERATION) 


Figure 23.- Longitudinal guidance laws 
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X 


ZC6S ■(X-XFLI) TAN (GS) + ZFLI 


ZCGS 



VELOCITY COMMAND LIMIT 


ZFLI -ALTITUDE AT INITIATION OF FLARE 
XFLI -RANGE AT INITIATION OF FLARE 


OUTPUT s 


DV2 WHEN INPUT >DV2 
INPUT WHEN 0.< INPUT <DV2 
0. WHEN INPUT <0. 


GS- GLIDE SLOPE ANGLE 

ZCGS -GLIDE SLOPE ALTITUDE AT CURRENT 
RANGE 


VELOCITY COMMAND RATE 

LIMIT *±DV5 ZFLI * - I50.FT 

NOTE : VELOCITY COMMAND LIMIT 

PROVIDES FOR A GLIDE SLOPE 
INTERCEPT 

DV2 * 15. FT/SEC 
DV 5 * 5. FT/SEC 2 


XFLI « - 1550. FT 
GS « I RAD 
DV4 * .2 SEC' 1 


Figure 24a.- 


Vertical guidance laws 


( | Z | > | ZFLI | ) 
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VERTICAL 

VELOCITY 

COMMAND 


vzco 



FLARE INITIATION 

TIME 


VZCO = 7.0 FT/SEC (INITIAL VELOCITY COMMAND FOR FLARE) 
TVZD = 23.0 SEC (TOTAL TIME OF OPEN LOOP FLARE) 


TAU IS ELAPSED TIME FROM FLARE INITIATION. WHEN TAU = TVZD 
THE SYSTEM REVERTS TO AN ALTITUDE HOLD SHOWN BELOW. 



VERTICAL 

VELOCITY 

COMMAND 


GVI = 0.1 SEC 1 
GV2 = 0.2 (FT/SEO/FT 


VZGR = (50.+ Z) GV2 (WHEN TAU s TVZD ) 


Figure 24b.- Vertical guidance laws ( f Z { < | ZFLI | ) 


SYSTEM SYNTHESIS 

The flight control modes were synthesized using classical 
linear analysis, fixed point hybrid simulation and mission 
simulation. First, gain estimates were made with linear analysis 
Then, the analysis was verified on the simulation. Based upon 
both the analysis and the simulation, final gain selection was 
made. This section presents the results of that synthesis 
procedure . 


Analysis 

This section presents the design development and system 
performance of each flight control mode, based on linear analysis 
Since the control system was to be implemented digitally, the 
effect of a sampled system as opposed to a continuous one had to 
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be considered. Based upon the recommendations of reference 9, 
a sample period of .032 sec. was initially chosen. This sample 
rate was considered high enough to warrant using continuous 
analysis in preliminary system synthesis and discrete system 
analysis for performance verification. This approach was taken 
with the digital system being modeled by a sample and hold network 
plus a computational lag in the forword loop. 

The design objective of each flight control mode was to fully 
develop its dynamic capability in realizing a desirable system 
transient response and to minimize the error of the output in 
following a step input. The step response specification was a 
fast response with as little overshoot as possible. Translated 
in terms of conventional performance criteria, the dominant 
system roots were desired to be slightly overdamped while making 
the system's bandwidth reasonably large with a minimum of compen- 
sation. Integral plus proportional compensation was inserted in 
the forward loops of several modes to assure zero steady state 
error . 

The CH-46C helicopter was modeled by two sets (longitudinal 
and lateral) of uncoupled differential equations with appropriate 
stability and control derivatives depending on the flight con- 
dition of the helicopter. The stability and control derivatives 
were obtained from the airframe manufacturer. The various flight 
conditions investigated are listed in Table I. The stability and 
control derivatives for Flight Condition 1 are listed in Appendix 
B. The perturbation equations of motion of the helicopter were 
written in terms of Euler angles and inertial velocities with 
respect to an Approach Navigation Frame (ANF) . These equations 
allowed the helicopter to be controlled with respect to a pre- 
selected landing site. In addition, each input to the airframe 
equations was passed through a second order actuator model with 
damping ratio 0.6 and frequency 15.0 rad/sec. The design approach 
for each flight control system was based on the normal flight con- 
dition (flight condition 1) at hover where the pitch axis is 
decoupled from the vertical axis and the roll axis is decoupled 
from the rudder axis. Compensation was first selected for this 
condition and then analyzed for each cruise flight condition 
where these aircraft axes became coupled. Adjustments were made 
in the few cases where the design was unacceptable at other 
flight conditions. 

The flight control modes are separated into the longitudinal 
axis and the lateral axis for convenience of presentation. 

Longitudinal Axis (General) 

The linearized longitudinal equations of motion of the 
helicopter in the ANF coordinate system are shown in Figure 25. 

The bare helicopter roots for flight condition 1 are shown in 
Figure 26 as a function of forward velocity. It is apparent 
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TABLE I.- CH46C AVAILABLE FLIGHT CONDITIONS 


Flight 
Condition # 

Gross Weight 
(LBS) 

Altitude 

(FT) 

C.G. 

Position 

Climb Rate 
(FPM) 

Forward Speed 
(KTS ) 

1 

13,400. 

Sea Level 

Normal 

0 . 

0,40,60,80,100,120,140 

2 



Most Fwd 



3 



Most Aft 



4 


10,000. 

Normal 


40,60,80,100 

5 



Most Aft 



6 

15,500. 

Sea Level 

Normal 

0 . 

0,40,60,80,100,120,140 

I 

7 



Most Aft 



8 

13,400. 

Sea Level 

Normal 

-1500. 

0,40,60,80 

9 

■ 


Normal 

+ 1500. 

.i 




that the helicopter exhibits an inherent instability throughout 
the flight velocity regime. The migration of the roots indicates 
a change from a pendulum type instability to a speed type in- 
stability (i.e., real root in right-half plane) as the helicopter 
passes from the hover phase through transition. 


[[(cm Os-f-wcot ®T« + 4r * ln fire)] [(-•'" Str-tH" 00 * O] [" 4? s + 4 ^ WtiTtJF u m + d co » ®t J Av x 4£ 

[(*'" ®ri») s ~(-ni Lco * ®m + Tjf*'" ®tr)] [( co# + » T „] AV Z * 4£ 8 

L[-T^ Ct »«TX-^* ,n «T„] [^* ine TX-^7 C0 * e Tx] [ SZ - S rS + S WT "'^ UT "] . A ®. X 

WHERE V %Tn & U TR cos 6^ + W TR sin 8 JK , B ■ I zz co* 8^- J„ (In 8 TR AND S* ^ THE LAPLACE DIFFERENTIAL 


OPERATOR. 


Figure 25.- Longitudinal perturbation equations 


80 

100 I 60 40 


- 2.0 


80 -- 7 

\ 40 * 

J00 


-ols 1 140 lo ,k) Ja5 

0 on 


Figure 26.- Longitudinal perturbation equation root migration 
as a function of airspeed (Flight Condition 1) 
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Attitude Command Mode I (Longitudinal Axis) 


The open-loop frequency response of the pitch to DEC channel 
of the helicopter at hover including actuator dynamics is shown 
in Figure 27. The poor low-frequency response prompted the 
insertion of an integral plus proportional compensation in the 
forward loop of the pitch axis. The integrator gain, GIG, is 
essentially the zero location of a lead term as shown by the 
block diagram of the integral plus proportional compensation in 
Figure 28. The optimum value of GIG was determined from the 
open-loop frequency response which is repeated in Figure 29. 

The results of making the break frequency, GIG of the lead term 
equal to the natural frequency of the open- loop system, .45 
radians/sec, are shown. The integrator increased the gain and, 
therefore, performance of the low frequencies while the lead 
term maintained marginal stability of the system. Making the 
break frequency greater than the natural frequency increased the 
bandwidth of the system but in the process destabilized the sys- 
tem. Making the break frequency less than the natural frequency 
reduced the destabilizing effect but the low frequency performance 
deteriorated. Selecting GIG = .2 radians/sec appeared to be a 
reasonable compromise. Note that the system is also gain in- 
sensitive as indicated by the fairly flat phase curve over the 
region of stability. With a lower break frequency on the lead 
term the system became somewhat more gain sensitive. 



- 60 1 


.03 


0.1 0.3 1.0 

FREQUENCY, RADIANS/SEC 


3.0 


7.0 


Figure 27.- Q/DEC open- loop frequency response at hover 
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PHASE, DEGREES 



MAGNITUDE, DB 


Figure 28.- Integral plus proportional compensation 



FREQUENCY, RADIANS/SEC. 


Figure 29.- 0/DEC open- loop frequency response at hover with 
integral plus proportional compensation in the 
forward loop as a function of GI0 


In order to achieve an adequately damped pitch attitude 
command system, some rate feedback was anticipated. A root- 
locus analysis was used to determine the rate feedback gain and 
attitude forward loop gain. The system configuration is shown 
in Figure 30. Figures 31, 32, and 33 show the attitude loop 


PHASE , DEGREES 




root locus with gain parameter G0 and with GQ = .25, .5, and 2. 

respectively. When GQ = .25 the locus emanating from the 
actuator poles was drawn into the real axis, but the locus 
emanating from the helicopter complex poles did not become 
sufficiently damped. Increasing GQ to .5 (this corresponds to 
moving the real open loop zero from -4 to -2) drew the complex 
locus to the real axis in the left half plane. This case was 
considered near optimum since the break-away angle from the 
complex poles is approximately 45° suggesting a maximum excursion 
of the locus into the left half plane. Increasing the rate gain 
from .5 to 2. drew the complex locus into the real axis abruptly 
resulting in a sluggish system. The attitude gain selected from 
Figure 32 was 20. inches/radian. 

Since the helicopter is not augmented in the vertical axis 
in Attitude Command Mode I, this completed the continuous analysis 
for this mode in the Longitudinal Axis. The system dynamics for 
Attitude Command Mode I were not verified analytically for the 
remaining flight conditions. This was left to the simulation. 

The root locus of the pitch attitude loop at hover with G0 
as the gain parameter and with sampler and hold plus a computa- 
tional lag in the forward loop is shown in the S plane in 
Figure 34. The system damping is reduced somewhat when the sys- 
tem is digitalized. 



Figure 30.- Pitch Attitude Command system 
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Figure 31b.- Root locus at origin of pitch loop at 
hover (G10 = .2, GQ = .25) 
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Figure 32a.- Root locus of pitch loop at 
hover (GI0 = .2, GQ = .5) 








Figure 34a.- Root locus of pitch loop at hover with sampler 
and hold plus computational lag in the forward 
loop (GIG = .2, GQ = .5). Sample period = .032 
second. 



Figure 34b.- Root locus at origin of pitch loop at hover with 
sampler and hold plus computational lag in the 
forward loop (GIG = .2, GQ = .5). Sample period 
= .032 second. 
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Attitude Command Mode II (Longitudinal Axis) 

The pitch axis for this mode is the same as in Attitude 
Command Mode I. A simple Rate of Descent Command system can be 
implemented in the vertical axis by feeding back the vertical 
velocity with forward loop gain as shown in Figure 35. The rate 
of descent Vz is defined positive downwards and collective rotor 
pitch (DCC) is defined positive for upward force or lift. There- 
fore, to achieve the correct sense in the simple rate of descent 
command system, the velocity gain GZD must necessarily be 
negative . 



Figure 35.- Simple rate of Descent Command system 


The root locus of the simple rate of descent loop with gain 
parameter GZD is shown in Figure 36. The output response speed 
can be increased as the velocity gain is increased. The closed 
loop frequency response corresponding to a velocity gain of 
-.2 in/ (ft/sec) is shown in Figure 37. The resulting bandwidth of 
the vertical axis is 3 radians/sec, while the steady state error 
was 20 percent. Increasing the velocity gain widened the band- 
width slightly and decreased the steady state error at the cost 
of decreasing the damping of the system. This simple system was 
thus precluded as an effective rate of descent command system. 

To improve the system, an altitude feedback loop was added 
in such a way as to be equivalent to inserting an integral plus 
proportional compensation. The system is shown in Figure 38. 

The root locus for this system with GZ as the gain parameter and 
with GZD = -.4 in/ (ft/sec) is shown in Figure 39. A reasonable 
choice for GZ , anticipating a decrease in damping with forward 
velocity, was 1.6 (ft/sec) /ft. The closed loop frequency response 
is shown in Figure 40. The steady rate of descent error was zero 
while the bandwidth was increased to 5. radians/sec. This 
represented a major improvement o TT er the previous rate of descent 
command system. Also, this system has altitude hold capability. 
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Figure 36a.- Root locus of the simple rate-of-descent loop 
at hover with pitch loop closed (GQ = .5, 

G0 = 20.) 
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MAGNITUDE, DB 


Figure 36b.- Root locus at origin of the simple 
rate-of-descent loop at hover with 
pitch loop closed (GIG = .2, GQ = 
.5, GG = 20. ) 



.03 0.1 0.3 1.0 3.0 7.0 


FREQUENCY, RADIANS/SEC 

Figure 37.- Vz/VzcoM closed-loop frequency response for simple 
rate-of-descent command system at hover with pitch 
loop closed (GIG = .2, GQ = .5, GQ = 20., GZD = -.2) 
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Figure 39b.- Root locus at origin of the altitude loop at hover 
with pitch loop closed (GI0 = .2, GQ = .5, 

G0 = 20. , GZD = -.4) 



Figure 40.- closed-loop frequency response at hover with 

altitude feedback and pitch loop closed (GI0 = .2, 

GQ = .5, G0 = 20., GZD = -.4, GZ = 1.6) 
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The root variations of the final Rate of Descent Command 
Mode system as a function of forward velocity are shown in 
Figure 41. The bandwidth and damping were both considered satis- 
factory for the flight velocities of interest. 
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Figure 41.- Root variations of the rate of descent Command 
System as a function of airspeed with pitch 
loop closed (GIG = .2, GQ = .5, G0 = 20., 

GZD = -.4, GZ = 1.6) 


The root locus of the Rate of Descent Command loop at hover 
with GZ as the gain parameter and with sampler and hold plus a 
computational lag in each axis of the Attitude Command Mode II 
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system is shown in the S plane in Figure 42. Only the damping of 
the system was changed (reduced slightly with GZ = 1.6 (ft/sec)/ft) 
as the system was digitalized. 



Figure 42a.- Root locus of rate-of-descent command loop 
at hover with sampler and hold plus compu- 
tational lag (GIG = .2, GQ = .5, G0 = 20. , 
GZD = -.4) Sample period = .032 sec 
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Figure 42b.- Root locus of origin of rate-of-descent command loop 
at hover with sampler and hold plus computational 
lag (GI0 = .2, GQ = .5, G0 = 20., GZD = -.4) Sample 
period = .032 sec 

Velocity Command Mode (Longitudinal Axis) 

The vertical axis of this mode is the same as in Attitude 
Command Mode II. A forward velocity command system was imple- 
mented by feeding back the forward velocity in the pitch axis as 
shown in Figure 43. Since positive pitch command generates 
negative V x the forward velocity gain must be negative. The 
root locus of the velocity command system at hover with velocity 
gain GVX as the gain parameter is shown in Figure 44. The 
velocity gain GVX = -.016 radians/ (ft/sec) was selected. The 
closed-loop frequency response is shown in Figure 45. It indicates 
a bandwidth of about 1 radian/sec. The velocity loop closure in 
the pitch axis had little effect on the vertical axis velocity 
response. The root variation of the Velocity Command Mode as a 
function of forward velocity is shown in Figure 46. The system 
was considered acceptable for all the flight conditions of 
interest . 

The root locus of the forward velocity loop at hover with 
GVX as the gain parameter and with sampler and hold plus a compu- 
tational lag in the forward loops of each axis of the Velocity 
Command Mode system is shown in the S plane in Figure 47. Again, 
the main effect of digitalizing the system is a slight change in 
the system's damping (GVX = -.016). 
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Figure 45.- Vx/VxC0M closed 
Velocity Coitunan< 
GQ = .5, G0 = 2 
GVX = -.016) 






Figure 46.- Root variation of Velocity Command System as a 
function of airspeed (GIG = .2, GQ = .5, 

G0 = 20., GZD = -.4, GZ = 1.6, GVX = -.016) 
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Figure 47.- Root locus of forward velocity loop of Velocity 
Command System at hover with sampler and hold 
plus computational lag (GIG = .2, GQ = .5, 

GQ = 20., GZD = -.4, GZ = 1.6) Sample period = 
.032 sec 
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Lateral Axis (General) 

The linearized lateral equations of motion of the helicopter 
in the ANF coordinate system are shown in Figure 48. The bare 
helicopter roots for Flight Condition 1 are shown in Figure 49 
as a function of forward velocity. The first task of the lateral 
system was to simply stabilize the vehicle. 


[ S -^-] [-£ S -(« C0 * «TR + TiT W Tr)] [(l£* in ®TR-^ C0 * ®Tr) S + T^'' xtr ] AV y ^ 

[ S ^ ~ S ~ ^TR~* in ^TR) "** ** n ^"R — ®Tr) ^ 1 

^W TR ] [s 2 ♦ (^.in e TR - £co. B rtt ) S ♦ £v xtr ] 

WHERE V XTR $ U TR co* fl TO + W TR sin fl TR . B $ cot « TR - J„*in « TR AND S - ^ THE LAPLACE DIFFERENTIAL 
OPERATOR. 


Figure 48.- Lateral-directional perturbation equations 



Figure 49.- Lateral-directional perturbation equation root 
migration as a function of airspeed (Flight 
Condition I) 


55 



Due to the near coincidence of the complex zeros with the 
complex divergent poles, as shown in Figure 50, it was apparent 
that the "Dutch Roll" mode could not be stabilized with simple 
yaw rate feedback through the rudder channel. However, the roll 
axis zeros shown in Figure 51 indicated that the divergent oscil- 
latory "Dutch Roll" mode could be stabilized and satisfactorily 
damped by applying roll attitude plus roll rate feedback through 
the DAC channel. The roll axis, therefore, was considered first. 



- 2.5 - 2.0 - 1.5 - 1.0 - 0.5 0 . 0.5 1.0 

Figure 50.- V/6 transfer function poles and zeros at 60 knots 



Figure 51.- $/6 a transfer function poles and zeros at 60 knots 
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The open- loop frequency response of the roll to DAC channel 
at hover including actuator dynamics is shown in Figure 52. The 
poor low-frequency performance indication prompted the insertion 
of an integral plus proportional compensation in the forward 
loop. This assured an ideal steady-state follower system. The 
open-loop frequency response as a function of integrator gain, 
GI$, is shown in Figure 53. GI$ was chosen to be .35 sec-1. The 
rationale for this selection is similar to that given in the 
longitudinal axis. 



FREQUENCY, RADIANS/SEC. 


Figure 52.- <£>/DAC open-loop frequency response at hover 


The roll attitude command system was then configured as 
shown in Figure 54. The roll rate gain TP was determined from 
the open-loop frequency response of the DAC to roll channel at 
hover with integral plus proportional compensation in the forward 
loop. This is shown in Figure 55. A value TP = .5 sec provided 
the potential for sufficient phase margin (system damping) without 
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degrading the system bandwidth. The roll angle gain, G$, was 
determined from a root-locus plot of DAC to roll channel with 
GI$ = .35 sec - -*- and TP = .5 sec and with G$ as the gain parameter. 
This is shown in Figure 56. G$ = 15 in/rad was selected as the 
roll angle gain. 



.03 0.1 0.3 1.0 3.0 7.0 

FREQUENCY, RADIANS/SEC 


Figure 55.- <J/DAC open-loop frequency response at hover with 
0.35 integral plus proportional compensation with 
varying rate feedback 


Note that with this gain value a slowly divergent heading 
mode exists and augmentation in the yaw channel was required. 
Stability was achieved with a yaw damper control loop as shown 
in Figure 57. The root locus of the rudder to yaw rate channel 
with GR as the gain parameter is shown in Figure 58. With 
GR = 30 in/ (rad/sec) the heading was stable and the "Dutch Roll" 
damping was not affected. However, examination of the mode 
during cruise revealed a root in the right-half plane. This 
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Figure 56a.- Root locus of roll attitude loop at hover 
(GI$ = 0.35, TP = 0.5) 



Figure 56b.- Root locus at origin of roll attitude loop 
at hover (GI $ = 0 . 35 , TP = 0 . 5 ) 
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instability was removed at high airspeeds with a sideslip feed- 
back loop in the yaw channel. The Sideslip Stabilization System 
during cruise is shown in Figure 59. The root locus of DRC to 
sideslip angle channel with G6 as the gain parameter is shown in 
Figure 60 for the 100-knot flight condition. G3 = 1. (rad/sec)/ 
rad was selected as the sideslip angle gain. Note that positive 
feedback was required. 




Figure 58a.- Root locus of yaw rate loop with roll loop closed 
at hover (GI$ = 0.35, TP = 0.5, G$ = 15.0) 
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Figure 58b.- Root locus at origin of yaw rate loop with roll 
loop closed at hover (GI$ = 0.35, TP = 0.5, 

G3> = 15.0) 



Figure 59.- Sideslip stabilization system during cruise 


Attitude Command Modes I and II (Lateral Axis) 

The mission requirement for the roll axis was roll attitude 
command. This was provided for in the basic stabilized system. 
The mission requirement for the yaw axis was for heading hold at 
low airspeeds and turn coordination at high airspeeds. To 
improve the turn coordination a cross-feed drive between the yaw 
rate and roll channels was developed assuming small angles and a 
simple roll angle response. The resulting cross feed-drive 
function became 


R = $ A $ 

C (TCFR)S + 1 C = (TCFR)S + 1 C 
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FIGURE 60b 


Figure 60a.- Root locus of sideslip feedback loop with roll and 
yaw rate loops closed at 100 knots (GI$ = 0.35, 


Figure 60b.- Root locus at origin of sideslip feedback loop with 
roll and yaw rate loops closed at 100 knots (GI$ = 
0.35. TP = 0.5. G$ = 15.0. GR 



GCFR and TCFR were evaluated at a nominal forward velocity of 
60 KTS giving GCFR = .3 and TCFR = .5 sec. The roll channel now 
controlled roll attitude while the crossfeed drive assisted the 
sideslip feedback in turning flight. The roll angle command 
system is shown in Figure 61. 



Figure 61.- Roll Attitude Command System for 
coordinated flight during cruise 


At low airspeeds, a heading hold command system was 
accomplished with the system depicted in Figure 62. The heading 
gain, Gip , was determined from a root locus of the heading channel 
at hover. The value Gip = 1.0 (rad/sec) /rad provided a reasonable 
trade off between the speed of the heading loop and the actuator 
damping . 


Velocity Command Mode (Lateral Axis) 

The velocity mode required a lateral velocity command 
system in the vertical heading axis system. This lateral velocity 
command mode was implemented with a velocity loop closure around 
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Figure 62.- Heading Command System at low speeds 


the roll axis with the yaw axis remaining as in the Attitude 
Command Modes. This is shown in Figure 63. The root locus of 
the lateral velocity command loop with GVY as the gain parameter 
is shown in Figure 64. GVY = .02 rad/ (ft/sec) was selected as 
the velocity gain. The root variations with airspeed of the 
Lateral Velocity Command System for Coordinated Flight, with 
GVY = .02 rad/ (ft/sec) , is shown in Figure 65. The roots of 
the system are adequately damped for all the flight velocities of 
interest . 



Figure 63.- Lateral Velocity Command Mode 


Automatic Mode (Lateral Axis) 

The lateral axis of the Automatic Mode was implemented by 
closing lateral position and velocity about the roll command sys- 
tem. This is shown in Figure 66. The root locus of this system 
with GY as the gain parameter is shown in the Figure 67. A gain 
value of GY =0.2 (ft/sec)/ft was selected. 
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Figure 64a.- Root locus of lateral velocity loop with roll, yaw 
rate, and sideslip loops closed (GI4> = 0.35, TP = 
0.5, G4> = 15.0, GR = 30., G3 = 1.0) at 100 knots 



Figure 64b.- Root locus at origin of lateral velocity loop with 
roll, yaw rate, and sideslip loops closed (GI$ = 
0.35, TP = 0.5, G4> = 15.0, GR = 30., G3 = 1.0) at 
100 knots 
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Figure 66.- Lateral position loop of the automatic modes 



Figure 67.- Root locus at origin of lateral position feedback 

loop with roll, yaw, sideslip, and lateral velocity 
loops closed at 100 knots 


Final System Definition 

With the linear analysis completed, the final system was 
configured and the final gain selections were determined from the 
simulation of the flight control modes. First, the agreement of 
the linear analysis and the simulation was confirmed via the 
correspondence of crossover gains and frequencies. Then, the 
performance of each mode was evaluated for the gains determined 
from the linear analysis. This evaluation was performed with 
both fixed-point and mission-simulation. Most of the gains were 
found acceptable with the exceptions noted below. 

In the longitudinal axis, the velocity system response was 
too fast for pilot acceptance (e.g., the helicopter pitched over 
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too far for small velocity commands) . The effective bandwidth 
of this system was therefore reduced by decreasing GZD, GZ , and 
GVX to -0.2 in/ (ft/sec) , 1.3 (ft/sec)/ft, and -.01 radians/(ft/ 
sec), respectively. In addition, an integral plus proportional 
compensation was found necessary in the forward velocity command 
system. This served as a trim adjustment on pitch attitude 
which is a nonlinear function of forward velocity. The integrator 
gain GVXI , was adjusted to -.001 radians/ (ft/sec) . 

In the lateral axis, the bandwidth of the yaw channel was 
too wide, allowing high-frequency wind noise to excite this axis 
excessively. To avoid this, GR and G3 were reduced to 15. 
in/ (radian/sec) and .3 (radian/sec) /radian respectively, and at 
high airspeeds a first-order lag (time constant T$ = .5 sec) was 
added in the sideslip feedback loop. On the other hand, the 
velocity command system was too slow, hence, GVY was increased 
to .04 radian/ (ft/sec) . 

The lateral axis of the automatic system was adjusted during 
simulation to provide maximum performance in the presence of 
gusts. Accordingly, at low airspeeds a crossfeed drive with 
large time constant was inserted to absorb steady state winds 
through a heading change while maintaining the localizer track. 
This system, together with the gain values, is presented in the 
block diagrams . 

Finally, all the stick sensitivities were determined 
experimentally and are listed in the block diagrams. 


SIMULATION DESCRIPTION 

The source of performance data for the various low-visibility 
landing systems was a fixed-base simulation. The simulation was 
comprised of a single-seat cockpit and a hybrid computer. The 
analog portion of the hybrid computer simulated the helicopter, 
while the digital portion simulated the digital multi-mode landing 
system. In this section, the physical characteristics of the 
simulator are presented; the simulated helicopter equations of 
motion are developed; the simulation set-up procedure is discussed; 
the gust model, which was implemented on the digital computer, 
is presented, and the technique used to determine the performance 
level of each approach is shown. 


Cockpit Description 

The simulation cockpit represented one-seat of a helicopter 
with the standard controls (i.e., pedals, center stick, and the 
collective stick to the left of the pilot) . The pedals and 
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center stick were force centered about their respective null 
positions. The force gradients were: 


Longitudinal Center Stick 
Lateral Center Stick 
Pedals 


Measured 

2 . 0 lb/ in 

1.0 lb/ in 

23.0 lb 


Acceptable Range (Ref. 10) 

0.5 to 2.0 lb/in 

0.5 to 2.0 lb/ in 

15.0 lb (max. deflection) 


Although the pedal forces were slightly high, they proved 
acceptable to the pilots. The center stick did not have a trim 
capability and the collective stick had an adjustable friction 
setting . 

The instrument panel layout (Figure 18) was based on the 
arrangement preferred by B.L.E.U., reference 11, with the addi- 
tion of both a range and a range rate indicator. The Flight 
Director Indicator (Figure 3) was used to display attitude 
information and as a three-axis "fly to" command display. An 
upward displacement of the vertical tab was interpreted as a 
command to fly up by means of the collective stick. A displace- 
ment to the left by the vertical bar meant that the cross-range 
rate was to be increased to the left by means of the lateral 
center stick. An upward displacement of the horizontal bar was 
interpreted as a command to pitch up by means of the longitudinal 
center stick. The radar altimeter was scaled from zero to 200 ft. 
The Mode Selection Panel consisted of two buttons which permitted 
the pilot to begin his approach after the completion of the 
simulation setup and to return the simulator back to the setup 
mode at the end of the mission. The digital range indicator had 
a resolution of .01 nautical miles, and its full scale was 9.99 
nautical miles. Figure 68 shows the overall cockpit. Although 
the instrument panel in the photograph is an earlier version of 
the one described, the two are similar. 


Analog Computer Hardware 

The cockpit interfaced with the Beckman 2200 analog computer 
which simulated the helicopter's dynamics, powered the cockpit 
instrumentation, and served as a transfer point for signals to 
and from the cockpit and the digital computer. The following 
list has been included to show the capability of the analog 
computer and the extent to which it was utilized. 
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Figure 68.- Simulator cockpit 
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Component 

Number Available 

Allocation 

Operational Amplifiers 

120 

100 

Integrators 

72 

36 & 29 as summers 

Multipliers 

48 

28 & 10 as summers 

Function Generators 

12 

10 

Bistable Multivibrators 

40 

10 

Monstable Multivibrators 

18 

5 

Or Gates 

114 

95 

Resolvers 

3 

3 


Approximately 80 percent of the analog computer's capacity, 
therefore, was utilized for this simulation. Of this percentage, 
60 percent was required to model the vehicle's dynamics and the 
remaining 40 percent was required to operate the cockpit. 
Associated hardware consisted of two eight-channel recorders, 
which were used to monitor pertinent state variables during each 
approach . 


Digital Computer Hardware 

The SDS 9300 digital computer interfaced with the analog 
computer and simulated the multimode digital flight control sys- 
tem, generated gusts, and reduced data pertinent to the performance 
of the helicopter during each approach. The computer cycle time 
was 1.75 p seconds, add time was 2 cycles for single precision 
and from 6 to 11 cycles for floating point precision, multiply 
time was 5 cycles for single precision and 8 cycles for floating 
point precision. The allocation of its memory was: 


Flight Control/Flight Director Laws 16K 
Guidance Laws 7K 
Performance Data 2K 
Gust Generation IK 
Miscellaneous 2K 
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This was a total of 28K or 90 percent of the computer's 32K 
storage capacity. This storage requirement could have been signif- 
icantly reduced if the program, which was written in a Fortran-like 
language, had been rewritten in assembly language. Associated 
hardware consisted of the typewriter, through which the operation 
of the simulator was controlled, and the line printer, for digital 
performance data output at the end of each approach. 


Interface 

The interface between the two computers consisted of analog 
to digital lines (ADL) , digital to analog lines (DAL) , interrupt 
lines, signal lines, and test lines. Interrupt lines permitted 
the analog computer to affect the digital computer by means of 
pulses, while the signal lines allowed the converse to take place, 
and the test lines permitted the digital computer to ascertain 
the state of certain components in the analog computer. The 
capacity of the interface and its utilization in this simulation 
are shown in the following list. 


Trunk Line 

Number Available 

Allocation 

ADL ' s 

40 

20 

DAL ' s 

40 

33 

Interrupt Lines 

10 

1 

Signal Lines 

16 

4 

Test Lines 

16 

11 


ANALOG COMPUTER PROGRAM 


Actuator System Equations 


The actuator model lumped the dynamics of the electro- 
hydraulic actuator, the lower hydraulic power boost, the upper 
hydraulic power boost and all the interconnecting linkage into 
a second order lag. The model was the same in each of the four 
input channels and the equations were: 
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( 1 ) 


A + 2T, ii) A + (jj A 
c a a c a c 


A + 25 to A + ( 1 ) A 
e a a e a e 


A + 2£ to A + to A 
a a a a a a 


A + 2? to A + to A 
r a a r a r 


to 2 DCC 
a 

2 

to DEC 

ci 

2 

to DAC 
a 

to 2 DRC 

a. 


where C a = 0.6 and to a = 15 radians/sec. Since each actuator 
travel is limited, the following limits were imposed on the out- 
puts from the actuator models. 


A 

c 

travel 

limited 

to 

0 

• 

to 

10.0 inches 

A 

e 

travel 

limited 

to 

± 

3 

.0 

inches 

A 

a 

travel 

limited 

to 

± 

3 

. 0 

inches 

A 

travel 

limited 

to 

± 

4 

.0 

inches 


Rotor System Equations 


Forces and moments on the helicopter are produced by the 
rotors. Either the collective angle of attack of the rotor 
blades or the cyclic angle of attack is changed by the control 
input. In either case, the rotor system exhibits some dynamics 
before settling out at the new operating point. In this simula- 
tion, the dynamics were modelled by a second order system. This 
model appeared in all four channels of input. 


The equations were : 


6 + ^ 6 + ft 2 6 

c 8 c < 


Q 2 A 


6 + ^ 6 + Q 2 6 

e 8 e ( 


2 

n a 


6 + ^ 6 + n 2 5 

a 8 a < 


ft 2 A 


«r + : fK + 


Q 2 A 


( 2 ) 
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where the Locke Number, y = 10. and the angular rate of the 
rotor ft = 27. radians/sec. 


Helicopter Equations 

The simulation equation for the forward velocity of the 
helicopter is developed in this section. The other five equations 
of motion were determined in a similar fashion. 

The nonlinear forward velocity equation is written: 

X 

U = RV - QW + — - g sin 0 (3) 

The terms in Eq. (3) were easily implemented with the exception 
of X^, which contains the aerodynamic characteristics of the 
helicopter. The aerodynamic force, X^, is a function of the 
helicopter state variables. It was assumed that the longitudinal 
and lateral-directional modes were aerodynamically decoupled 
and, therefore, that: 


X A = X A (U, Q, W, 6 e , 6 c ) . (4) 

Although not always a valid assumption for helicopters, it was 
believed to be valid for the CH-46C. This assumption had been 
made with success by other investigators for this helicopter. 

The force was then expanded about a steady flight condition 
of level flight. For such a flight condition, the specification 
of forward velocity specified Q TR = O, W TR (U TR ) , <$eTR (^Tr) , 
6 c tr( u Tr)' an< 3 ®Tr( u TR^* The expanded force was then: 


X A (U 

TR ,0,W TR^ U TR^ ' 6 eTR * U TR* ' 6 cTR (U TR ) ) 


3X. 



A 

9U 

^ U TR'° ,W TR (°TR^ ' 6 eTR ^ U TR^ ' 5 cTR (U TR^ ^ 

(U-U TR ) + 

3X, 



A 

96 

^ U TR ,0,W TR ^ U TR^ ' 5 eTR ^ U TR^ ' 5 cTR ^ U TR^ ^ 

^ <S c“ lS cTR ) 


In using the expanded force in the simulation, the trim point was 
varied as the forward speed of the aircraft varied. In this way. 
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the forces would always be expanded about a steady level flight 
condition whose velocity would be the present velocity of the 
helicopter. Thus, in Eq. (5) = U dropping the second term 

to give: 

X A =■ X A (D,0,W TR (D),« eTR (0),6 cTR (U)) 


3X, 


Tjw < u '°' w tr (u) 'W u, ' 6 ctr < u)> < k - w tr ( u)) 


( 6 ) 


g X 

+ nr <0 ' O ' W TR (U) ' 6 eTR (a) ' S C TR (U) > < V^TR < U > > 


The partials are the stability and control derivatives as a 
function of forward flight speed and were directly available. 

To construct the first term, the total derivative was taken with 
respect to U to give: 


dX 


A 


dU 


(U,0,W TR (U) ,6 

ax. 


eTR (U) ' cTR 


8 X 

~ | (U,0,W rpD (U),6^ D (U), 


TR 


eTR 


6 cTR (U ^ + 9W m „ (U, ° ,W TR (U) ' 6 eTR (U) ' 6 cTR (U) ) — ™ (U) + ••• 


TR 


9U 


9x a 86 ctr 

96 cTR (U,0,W TR ' 6 eTR (U ^ ' 6 cTR ~9U 


The lead term is simple and was an available stability derivative, 
while the partials with respect to the variables W TR (U), 6 c Tr( u ) 
6 s tr(U) were not available. Consequently, as a first approximation 
of X^ it was assumed that trim variations with forward speed were 
slow and that the products of partials in Eq . (7) were small when 

compared to the lead term. The first term of Eq . (6) was then 

constructed by integrating the X u stability derivative from 0 
to U to give: 

/ d 3 y 

W (o»0,w <a), 

' J ( 8 ) 


W° ) 'W 0 > ) d ° + X A (0 '° 0 - 4 eTR (0) ' «cTR (0) > 
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The wind gusts, Gxa and G ZA are a l° n 9 the body axes, X B and 
Z B . 


In principal, the forward velocity equation could now be 
implemented on the analog computer. Hardware limitations, however, 
required further equation modification. In particular, the 
number of function generators available to simulate the velocity 
dependence of all the stability and control derivatives and trim 
variables was not sufficient. As a result, an analytic study 
was made to determine which derivative changes had only a minor 
influence on the helicopter's aerodynamic characteristics as a 
function of forward velocity. It was found that twenty of the 
thirty derivatives could be satisfactorily held at their hover 
values. These derivatives are listed in Table II. Figures 70(a) 
and (b) present the plots of the remaining velocity dependent 
derivatives . 


Zw/m Zgc /m Mw/Iyy 



AIR SPEED AIR SPEED AIR SPEED 

kts kts kt$ 


M 8e /I yy 



AIR SPEED 
kts 


Figure 70a.- Air speed dependent coefficients 
used in the simulated helicopter 
equations of motion 
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TABLE II.- VALUES OF THE AIR SPEED INDEPENDENT COEFFICIENTS 
USED IN THE SIMULATED HELICOPTER EQUATIONS OF 
MOTION 


V" 

.05449 (f t/sec^) / (ft/sec) 

X q /m 

.60185 (ft/sec 2 )/ (rad/sec) 

X 6e /m 

.17696 (ft/sec 2 )/in 

X 6c /m 

1.2048 (f t/sec 2 )/in 

Z q/ m 

-.71511 (ft/sec 2 )/ (rad/sec) 

Z 6e /m 

-.00407 (ft/sec 2 )/in 

M /I 

q yy 

-.73173 (rad/sec 2 )/ (rad/sec) 

/I 

oc yy 

-.04765 (rad/sec 2 ) /in 

Y v /m 

-.02664 (ft/sec 2 )/ (ft/sec) 

Y p /m 

-.76514 (ft/sec 2 )/ (rad/sec) 

Y r /m 

-.12517 (ft/sec 2 ) / (rad/sec) 

W” 

.99794 (ft/sec 2 )/in 

Y 6r /m 

.14652 (f t/sec 2 ) /in 


-.50730 (rad/sec 2 )/ (rad/sec) 

V J X* 

-.02297 (rad/sec 2 )/ (rad/sec) 

L 5a /,I xx 

.46536 (rad/sec 2 )/in 

L. /I 
or xx 

-.12638 (rad/sec 2 )/in 

N /I 
p' zz 

-.01831 (rad/sec 2 ) / (rad/sec) 

N «a /r zx 

.03001 (rad/sec 2 ) /in 


.17584 <rad/sec 2 )/in 
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Lv/Ixx 



AIR SPEED 
kts 


Nv/X zz Nr/I zz 




AIR SPEED 
kts 


AIR SPEED 
kts 


Figure 70b.- Air speed dependent coefficients used in 

the simulated helicopter equations of motion 


As one further step in reducing the number of function generators 
required, the trim variables W TR (U), S c tr(U), ^eTR^ were also 
set to their hover values. The effects of neglecting these 
variations is discussed later in the section. The final form 
of the forward velocity equation used in the simulation is pre- 
sented along with the other helicopter equations in Eq . 10. 


“ = k [X A (U AS )+X q (0)Q+X w <0)(W AS- W TR (0))+X 6e (0)(6 e- 6 eTR <0)) 
+ X 6c (0) (<S c - 6 cTR (O) ) ] + RV - QW - g sin 0 


* - E [ V D AS )+ V° >Q+ V U AS> < W AE- W TR (0) > +Z 6e (0) ( V { eTR (0 > > 


+ Z 6c (0) (6 c -6 cTR (°) ) ] + QU - PV + g cos $ cos 0 
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Q - (M A <U AS )+M q< 0 > Q+M K <U AS )(W AS- W TR <0)) 

+ M 6e (U AS> <5 e- 5 eTR (0),+M 6c (0) <' “c'^cTR <°> > 1 


v “ 5 [ Y v (0)V AS +Y P (0)p+Y r <0)R+Y 6a< 0, ' 5 a +Y 6r (0)5 r I 


+ PW - RU + g sin $ cos 0 


p - TIT [I xz R+ V 0 )p+L r< o > R+ V u AS> v AS 

XX 


+ L, (O) 6 +L„ (O) 6 ] 

6a a 6r r 


R = 


'ZZ 


[ I xz P +Np (O) P + N r (U AS ) R+N v (U AS ) V AS 


+ N x (O) 6 +N X (O) 6 ] 
6a a 6r r 


( 10 ) 


The helicopter variables along the moving trim are defined as : 


U TR = U AS 


S eTR (U AS> 


W°> 


V TR (U AS ) ° 


<S cTR (U AS ) 6 cTR (0) 


W TR (U AS } = W TR (0) = ° 


6 aTR * U AS ) 


O 


P TR (U AS ) 


O 


6 rTR (U AS^ ° 


q tr (u as ) 0 


r tr* u as^ 0 


( 11 ) 


The wind gusts along the X B , Y B , and Z B body axes are G x ^, Gy^, 
and G za , respectively. The aircraft to air mass velocities are 
defined as 
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U AS 

= U 

g xa 

V AS 

= V 

- g ya 

W AS 

= w 

g za 


( 12 ) 

The state of the helicopter at the beginning of each approach is 


g 

H 

o 

II 

G 

O 

P IC 

w 

IC 

O 

s 

II 

Q IC 

< 

H 

o 

= 0 

R IC 


(13) 


The approximations made during the development of the 
simulation version of the helicopter equations were verified. 

Two approximations that were made during the equation develop- 
ment would result in an unknown degree of simulation trim error. 
First, the trim terms X A , Z A , and M A were constructed by in- 
tegrating the X u , Z u , and M u stability derivatives respectively 
from 0 to U. The second approximation came about when the 
velocity dependent trim variables W TR (U T p), <S e TR ' ^cTr( u tr) 

were set to their hover values. To check the net effect of these 
two approximations on trim, simulation trim values were generated 
and compared against their true values. Figure 71 shows the 
comparison. Agreement between the simulated and the actual trim 
values was satisfactory. 

Next, the approximation of the majority of the velocity 
dependent stability and control derivatives by their hover values 
would have an effect on the simulated aerodynamic characteristics 
of the helicopter. The primary aerodynamics of the uncontrolled 
vehicle are determined by the pole-zero locations of the open 
loop transfer function associated with each input channel, namely: 
0/6 e , W/6 C , $/<5 a , and R/6 r . Figures 7 2 through 77 show a comparison 
of the aerodynamics (location of the poles and zeros) of the 
simplified aircraft model used in the simulation with those of the 
original model for each of these transfer functions. There was 
satisfactory agreement with respect to both the actual location 
of the roots and their migratory trends as a function of forward 
velocity. The exception to this general agreement was in the case 
of the numerator associated with W/S c (Figure 74) . Although the 
simplified model obviously did not conform to the original in 
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Figure 71.- Comparison of simulation trim values 
with the actual trim values for the 
CH-46C helicopter 
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this case, the general overall agreement of the simplified model 
with the original led to the acceptance of this discrepancy. 
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Figure 7 2.- The roots of the denominator associated with 
the longitudinal motions of the CH-46C 
helicopter as a function of air speed in 
knots 
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OPERATING CONDITION 

GROSS WEIGHT : 13400 LB (NORMAL ) c.g POSITION : NORMAL 

ALTITUDE SEA LEVEL CLIMB RATE : 0.0 fpm 


Figure 73.- The roots of the numerator associated with the 

pitch channel (attitude) of the CH-46C helicopter 
as a function of air speed in knots 
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Figure 74.- The roots of the numerator associated with the 
power channel (vertical velocity) of the CH-46C 
helicopter as a function of air speed in knots 
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Figure 75.- The roots of the denominator associated with the 
lateral motions of the CH-46C helicopter as a 
function of air speed in knots 
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GROSS WEIGHT: 13400 LB (NORMAL) c.g. POSITION : NORMAL 

ALTITUDE: SEA LEVEL CLIMB RATE: 0.0 fpm 


Figure 76.- The roots of the numerator associated with the 

roll channel (attitude) of the CH-46C helicopter 
as a function of air speed in knots 
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Figure 77.- The roots of the numerator associated with the 
yaw channel (yaw rate) of the CH-46C helicopter 
as a function of air speed in knots 
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Euler Angle Equation s.- The following equations used body 
rates to generate Euler angles. 

0 = Q cos $ - R sin $ 


* 


sin $ _ cos $ 

cos 0 cos 0 


i = P + 'i' sin 0 


(14) 


The value of these variables at the beginning of each approach 
were 


IC 


0 TR (U IC ) 




IC 


*0 


$ 


IC 


o 


(15) 


Vertical H eading (VH) — Coordin ate System V e locity 
Equations . - The following equations resolved the body referenced 
velocities through 4> and 0 to give what was called VH-coordinate 
velocities . 

V^ = U cos 0 + V sin $ sin 0 + W cos $ sin 0 

A 

Vy = V cos $ - W sin $ 

v_ 

VI 1 = -u sin 0 + V sin $ cos 0 + W cos $ cos 0 

(16) 

Approach Navigation Frame (ANF) — Coordinate Sys te m Ve 1 o c i ty 
Equations . - The following equations resolved the VH-coordinate 
velocities through 'H to give the ANF — coordinate velocities. 

V x = V* cos Y - Vy sin Y 

V v = v5 sin Y + V^ cos ¥ 

i A Y 
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DIGITAL COMPUTER PROGRAM 

The digital program had the tasks of setting up the simulator 
for each run, simulating the digital flight control system, 
generating gusts, and reducing the data generated during each 
approach . 

The simulation set up prior to each approach was carried 
out by means of the typewriter I/O device. The setup consisted 
of specifying: 

(1) The flight control mode that was to be flown during the 
next approach (e.g.. Attitude II without Flight 
Director Commands) . 

(2) The state of the helicopter at which the approach was 

to begin (i.e., X IC , Y IC , Z IC , (vh) IC , (vh) IC , and y IC ) . 
The rates. Pic, QiC' an< ^ R IC ' were zero at approach 
initiation . 

(3) The direction of the mean wind vector in the ANF 
coordinate system. 

(4) Any changes with respect to the nominal control system 
gains, trajectory characteristics, or wind character- 
istics. This option was particularly convenient during 
the verification of the linear analysis and in 
establishing the final form of the flight control sys- 
tem and of the mission. 

At this point, the helicopter was commanded to fly to the 
specified initial conditions with a modified automatic control 
system. This procedure automatically established the initial 
conditions on the associated integrators and function generators 
in the analog computer. After the helicopter reached the starting 
condition, the pilot would begin the approach with the specified 
control system. At the completion of the mission, the simulator 
entered the set-up mode once again. 

In order to simulate a sample data flight control system, 
the digital program was structured such that it consisted of an 
idle loop and one major subroutine in which the flight control 
computations were carried out. Upon completion of this sub- 
routine, the program remained in the idle loop until a timing 
pulse was received from the clock in the analog computer. The 
interval between pulses (32 msec) represented the sampling 
period. 

In order to exercise the landing system modes, it was 
necessary to simulate winds. The model chosen was simple but 
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sufficient for the purpose of the experiment. Model parameters 
were taken from reference 12. The wind axis system was defined 
with the X axis along the mean wind vector, the Y axis in the 
plane of the local horizontal, and the Z axis down the local 
vertical. The wind model consisted of a mean-wind vector plus 
gusts along the wind axis. The gusts were assumed uncorrelated 
and consisted of zero-mean white Gaussian noise passed through 
a simple first-order filter. A mean-wind amplitude of 20 fps 
was chosen. This led to standard deviations for the X, Y, and Z 
wind axes gusts of 2.3, 1.6, and 1.0 fps, respectively. A filter 
time constant of 1.5 sec was chosen for each axis. 


In the simulation, the direction of the mean-wind vector 
could be varied from run to run while its amplitude was kept 
fixed. The total winds were then resolved from wind axes to 
body axes. The winds were added to the inertial velocities and 
passed through the stability derivatives to generate the forces 
and moments. A sample of the standard winds described above is 
shown in Figure 78. The derivation of the noise model is given 
in Appendix C. 


The digital computer reduced the performance data from each 
run. To establish system performance, the mission was divided 
into four phases, Localizer Track, ILS Beam Track, Flare, and 
Terminal Approach. During each phase, the digital computer 
recorded data consisting of the initial state, the absolute 
maximum deviation of certain variables from their nominal values, 
and, in certain cases, the root-mean-square values of certain 
deviations. At the end of each approach, a normalized performance 
index involving seventeen of these data points was calculated and 
printed out along with the data associated with each phase. A 
typical printout is shown in Figure 79. 

The form of the performance index was 


P . I . 


1 

17 


V (E ( 1 ) -ENQ M(i) ) HAX 
A, - EMAX (i) 


(18) 


where E represents each variable; ENOM represents the nominal 
value of each variable and EMAX represents the maximum allowable 
deviation of each variable from its nominal. Each of these terms 
was used in two ways to determine system performance. To establish 
the basic acceptability of each of the candidate control systems 
each of the terms associated with maximum absolute excursions 
from nominal (i.e., excluding the root-mean-squared terms) were 
looked at individually. If any term exceeded unity, the system 
was considered to have failed to meet the minimum performance 
specifications during that approach. The overall percentage of 
failed approaches to the total number of approaches made with a 
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Figure 78.- Example of the wind used in the simulation 
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MULTI-MODE HELICOPTER FLIGHT CONTROL SYSTEM 


INITIAL 

STATE 

MAXIMUM VALUES 


LOCALIZER TRACK _ 

X = -0 . 1 0000E 05 
Y= 0.00000E 00 

AY= 0.61035E-02 
AZ=-0 . 32080E 02 

ABYS= 0.89879E 00 
ABSDZ= 0. 12609E 01 

ABSAZ= 0.35132E 02 
ABSPHI= 0 . 1 9384E-01 


Z = -0 . 80000E 03 
VL0CE= 0.00000E 00 

PHI = 0.85205E-03 
THET A= 0. 49632E-01 

ABSVL0CE= 0 . 1 1850E-03 
ABSDVX= 0.13387E 01 

THET AP0S= 0.521 88E -01 
THETANEG= 0.49419E-01 


VX = 0.70982E 02 
VY=-0.31373E 00 

PSI=-0. 16508E 00 
BET A=-0 . 29996E-02 

ABSVY= 0 . 3704 1 £ 00 
ABSVZ= 0.79468E 01 

A B S P S I = 0.21599E 00 
ABSBET A= 0.1 4227E 00 

TT= 0.56608E 02 

VZ = 0 . 48828E-0 1 
AS= 0.87495E 02 

P=-0.21301E-03 
Q=-0 . 42603E-3 

ABSDAS= 0.12634E 02 
- AX= 0.16479E 01 

ABSP= 0.451 59E-01 
ABSQ= 0 . 29609E-01 

PLOCE= 0 . 41 21 OE-04 
PVEL= 0.48783E 00 

AX = -0 . 1 4771 E 01 

R= 0. 1 0651 E-02 

ABSAY= 0.26062E 01 ABSR= 0.38555E-01 

ILS BEAM TRACK 

PAS= 0.43343E 01 

X= -0 . 59981 E 04 
Y = - 0 .21 082E 00 

AY= 0.24414E 00 
AZ= -0 . 3 1 024E 02 

ABSY= 0.66576E 00 
ABSVL0CE= 0.17733E-03 

ABSAZ= 0.35608E 02 
ABS PH I = 0, 22579E-01 


Z=-0 , 59739E 03 
VL0CE= o . 30663E-04 

PH I = 0.5751 3E-02 
THETA= 0 . 1 3633E-01 

ABSGSE= 0.55401E-03 
A8SDVX= 0 . 97023 E 00 

THETAP0S= 0 . 39620E-01 
THETANEG= 0 . 1 256RF-01 

TT= 0.62464E 02 

GSE = 0 . 98547E-04 
VX = 0.70592E 02 

PSI=-0.17531E 00 
BETA=-0 . 51 41 4E -01 

ABSVY= 0.42948E 00 
ABSVZ= 0.90088E 01 

ABSPH I = 0.20407E 00 
ABSBETA= 0 . 97461 E-01 

PL0CE= 0.121 57E-03 

VY= 0.11985E 00 
VZ = 0.73120E 01 

P= 0 . 97986E-02 
0= 0 . 29609E-01 

ABSOAS= 0.11487E 02 

-AX= 0.13184E 01 

ABS P= 0. 34082E-01 
ABSO= 0 . 30248E-01 

PGSE= 0.341 78E-03 
PVEL= 0.81904E 00 

AS = 0.95564E 02 
AX= 0.21973E 00 

R= 0 . 29822E -02 

ABSAY= 0.18738E 01 

ABS R= 0 . 291 83E-01 

PAS= 0.10746E 02 

X = -0 . 1 5522E 04 

AY= 0.67139E 00 

FLARE 

ABSY= 0 . 2591 OE 01 ABSPHI= 0.38981 E-01 


Y = 0.61 465E-02 
Z = -0 . 1 4998E 03 

AZ=-0 . 31 354E 02 
PHI=-0.46864E-02 

Z = -0 . 51 91 3E 02 
ABSVLOCE= 0.1621 4E-02 

THETAPOS= 0.15848E 00 
THETANEG=-0 . T 81 06E-01 


VL0CE=-0.19927E-05 
GSE = -0 . 47340E-03 

THETA= 0.181 06E-01 
PSI =-0 . 1 751 OE 00 

ABSVY= 0.54204E 00 
ABSVZ= 0.78125E 01 

ABSPS I = 0.39088E 00 
ABSBETA= 0.17584E 00 

TT= 0.43200E 02 

VX = 0.75680E 02 
VY = 0. 33005E-01 

BETA= 0.181 60E-01 
P= 0. 63904E -02 

-AX= 0.43579E 01 
ABSAY= 0.14465E 01 

ABSP= 0. 35786E-01 
ABSQ= 0 . 41 964E-01 


VZ = 0.76782E 01 
AS = 0.86042E 02 

q= 0 . 48993E -02 
R= 0 . 1 0225E-01 

ABSAZ= 0.35675E 02 

ABSR= 0. 32378E-01 

PLOCE= 0. 1 2592E-02 

AX=-0. 40283E 00 


TERMINAL APPROACH 

[RANGE=1 00 . FT] 


X=-0.99742E 02 
Y« 0.13989E 00 

P H I = - 0 . 74554E-02 
THETA= 0. 1 26 1 OE 00 




Z = -0 . 51 894E 02 
VX= 0.15445E 02 

PS I = -0 . 38981 E 00 
BETA=-0 . 44465E -01 




V Y= 0 . 52847E-02 
VZ = 0.61035E 00 

P= 0. I 0438E-01 
Q=-0 . 63904E -03 




AS = 0.31159E 02 
AX = -0 .41 382E 01 

R= 0 . 21 301 E -02 
BEARING= 0 . 342 I 5E-03 




A Y = - 0 . 6 1 035E-02 
AZ=-0 . 32898E 02 

BEARING DESI RED=-0 . 1 3089E-02 
WORKLOAD TIME = O.OOOOOE 00 



CRITERIA 

LOCALIZER TRACK ILS MEAN TRACK 

FUNCTION 

FLARE 

TERMINAL APPROACH 

ABS DZMAX/ 1 00 . = 0.1 2609E-01 ABSGSEMAX/ . 035= 0.15829E-01 

ABSLOCEMAX/ . 1 = 0.11850E-02 ABSLOCEMAX/ . 1 = 0.17733E-02 

ABS [ZMIN + 50 . ] / 2 0 . = 0.95664E 
ABSLOCEMAX/. 1= 0.16214E-01 

■LOl ABSY/1 00 .= 0 . 1 3089E-02 

ABS fZ+50. 1/20.= 0. 94688E-01 

PLOCE/ . 05= 0 . 8241 9E-03 PL0CE/.05 

P VEL/20 . = 0 . 24391 E-01 PGSE/.018 

= 0.24313E-2 

= 0. 1 8988E-01 

PLOCE/. 05= 0.251 84E-01 

ABS [VX-1 6 . ] / 1 6 . = 0 . 34664E-01 
ABSBEARE/ [PI/61 = 0.31 549E-02 • 


PVEL/20 . = 

0 . 40952E-01 


ABSVZ/3 . = 0.20345E 00 


[SUMMATION OF THE TERMS / 1 7 

.= 0.34901-01 


Figure 

79.- Example 

of the digital 

performance data 

printed out at 


the termination of each simulated approach 



particular flight control system determined its acceptibility 
as an operational system. The Performance Index (P.I.) is a 
qualitative measure of over-all system performance which was used 
to rank the seven candidate control systems with respect to 
performance. It is not to be considered, in any sense, an 
absolute scale. 

The details of the P.I. can best be understood in the context 
of the nominal simulated mission. The division between the four 
phases of the mission was based on basic shifts in the set of 
guidance parameters that the pilot controls during the approach. 
During each phase, the associated terms in the P.I., when possible, 
consisted of these same guidance parameters, Table III. During 
the following description of how the P.I. relates to the mission, 
reference may be made to Figure 1, which is a diagram of the 
mission, and to the section entitled Low-Visibility Mission, in 
which the simulated mission is established and discussed. The 
Localizer Track phase was initiated at a range of 10,000 ft, an 
altitude of 800 ft, a ground speed of 42 kts (70 fps) , zero rate 
of descent, headed into the mean wind, aligned with the localizer, 
and approaching the glide slope intercept point. The pilot 
initially was concerned with maintaining altitude, ground speed, 
and localizer track. Consequently, for this phase to be con- 
sidered acceptable, the P.I. required that the altitude be 
maintained between 700 and 900 ft, that the ground speed root-mean- 
square deviation from 70 fps be less than 20 fps, and that the 
localizer deviation be within ± 0.1 rad (i.e., full-scale localizer 
sensitivity) , and that the associated localizer root-mean-square 
deviation be less than 0.05 rad (i.e., half-scale localizer 
sensitivity) . Somewhere between the range of 9000 and 8000 ft, 
the pilot would normally switch from his altitude hold task and 
commence to capture the glide slope. For this reason, the 
calculation of the associated altitude term in the P.I. was 
terminated at 9000 ft of range. In addition to assuming that 
glide slope capture would be initiated at some range less than 
9000 ft, it was further assumed that glide slope capture would be 
completed by the time the helicopter was at 6000 ft of range. 

Since the helicopter can acceptably approach the glide slope at 
an altitude of 700 ft, 1000 ft of range beginning at the 7000-ft 
glide slope intercept range was allowed for the intercept maneuver. 

At 6000 ft, the calculation of those terms associated with 
the Localizer Track were replaced by those of the ILS Beam Track. 
During this phase, the pilot would maintain the localizer, glide 
slope, and ground speed. The P.I. terms again monitored the 
deviations of these three parameters. The glide slope deviation 
from 0.1 radian was to be kept within ± .035 rad (2 deg) which 
again was the full-scale glide slope sensitivity. As the 
helicopter approached the flare/deceleration altitude of 150 ft, 
the pilot would progressively transfer his attention from main- 
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TABLE III.- GUIDANCE PARAMETERS 


Phase of Mission 

Parameter 

Nominal Value 

Perf. Index Parameter 

Localizer Track 
(From range 10,000 ft 
to range 6,000 ft) 

Altitude 

(terminated at range 
9,000 ft) 

800 ft 

Max . Error 


Localizer 

0. rad 

Max. Error 


Ground Speed 

70. fps 

RMS Error 
RMS Error 


ILS Beam Track 
(From range 6,000 ft 
to alt. 150 ft) 


Glide Slope 

(Term, at alt. 200 ft) .1 rad 


Localizer | 0. rad 

Ground Speed 70. fps 


Max . Error 
RMS Error 

Max. Error 
RMS Error 
RMS Error 


± 100. ft 


.1 rad 
.05 rad 
20. fps 


.035 rad 
.018 rad 

.1 rad 
.05 rad 
20. fps 


‘Tare 

(From alt. 150 ft 
to range 100 ft) 


Altitude 

Localizer 


50. ft 

(after flare) 
0. rad 


Lowest Pt 


Max . Error 
RMS Error 


± 20. ft 


.1 rad 
.05 rad 


Terminal Approach 
(From range 100 ft 
to touchdown) 


Altitude 

(at range 100 ft) 

Lateral Pos. 

(at range 100 ft) 

Course 

(at range 100 ft) 

Ground Speed 
(at range 100 ft) 

| Vertical Vel. 

(at range 100 ft) 


50. ft 
O.ft 


Aligned with center Error 
of pad 


16 . fps 
0. fps 


± 20. ft 

± 100. ft 











taining the glide slope to monitoring the radar altimeter in 
order to coordinate his stick responses with altitude. Con- 
sequently, the calculation of the glide slope terms was terminated 
at 200 ft of altitude. The Flare Phase was entered at 150 ft at 
which the pilot commenced to decelerate to a ground speed of 10 
kts (16 fps) and to flare to an altitude of 50 ft while maintaining 
the localizer track. The pilot then maintained a ground speed of 
10 kts, an altitude of 50 ft, and localizer track until a range 
of 100 ft. During the flare to the 50-ft altitude and subsequent 
altitude hold, safety considerations dictated that the helicopter 
should never drop below 30 ft. This established a lower altitude 
limit on flare. To establish an upper limit, a symmetric, 
although minimal, requirement would be that the helicopter must 
drop below 70 ft at some point. The associated altitude term in 
the P.I. represented these two conditions. At the range of 100 
ft, the Flare Phase terminated and the Terminal Approach phase 
was entered. At this range, it was assumed that the pilot would 
be able to switch to VFR for the final approach and touchdown, 
even in the worst weather conditions. Consequently, the simulated 
approach was terminated at this range; and the associated terms 
in the P.I. attempted to establish if the state of the helicopter 
at flare termination (i.e., breakout) would permit the pilot a 
high probability of completing the final approach and touchdown 
smoothly, without the need for abrupt maneuvers. For breakout 
to be considered successful, the helicopter's terminal altitude 
had to be between 30 and 70 ft, lateral position had to be less 
than ± 100 ft, course had to be within ± 30 deg of being aligned 
with the center of the pad, and the ground speed had to be less 
than 20 kts. One additional constraint was that the vertical 
velocity be less than 3 fps. This may seem overly restrictive, 
but it was found that whenever a flight control system repeatedly 
failed this test, the pilot's collective stick movement was 
excessive and that his ability to maintain a constant vertical 
velocity and altitude was unacceptable. 


EXPERIMENT 

The objective of the experiment was to determine the 
relative performance of the various landing modes. To do this 
six subjects flew a standard sequence of landing missions while 
the performance was recorded. The subjects consisted of four 
instrument-qualified pilots and two engineers with simulator 
experience. A summary profile on each of the six subjects used 
in the experiment is presented in Table IV. Subjects 1 and 2 
were engineers with no flight experience. They both had worked 
on the system development and understood the mission thoroughly. 
Each test sequence was run in two consecutive days. The first 
day was considered training and qualification on each mode. On 
the second day the production runs were made. 
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TABLE IV 


BACKGROUND OF EACH TEST SUBJECT 


IFR Certification 

1 

2 

3 

4 

5 

6 

Fixed When Last Certified 

Never 






Wing Duration of Certification 







Rotary When Last Certified 

Never 





' V 1 

Wing Duration of Certification 






■ 

IFR Instructor Certification 







Fixed When Last Certified 

Never 

Never 

Current 

Current 

Current 

Never 

Wing Duration of Certification 



7 yrs 

3 yrs 

2 yrs 


Rotary When Last Certified 

Never 

Never 

Never 

Current 

Never 

Never 

Wing Duration of Certification 

■ 



3 yrs 



VFR Instructor Certification 







Fixed When Last Certified 

Never 

Never 

Current 

Current 

Current 

Never 

1 Wing Duration of Certification 

1 



7 yrs 

6 yrs 

4 yrs 


1 Rotary When Last Certified 

■ Never 

Never 

Current 

Current 

Never 

Never 

j Wing Duration of Certification 



4 yrs 

6 yrs 



IFR Experience 







Fixed Wing 

0 hrs 

, 0 hrs 

| 150 hrs 

200 hrs 

200 hrs 

200 hrs 

Rotary Wing 

0 hrs ; 0 hrs 

] 0 hrs 

150 hrs 

0 hrs 

200 hrs 

VFR Experience 



1 




Fixed Wing 

0 hrs 

0 hrs 

1300 hrs 1000 hrs 1300 hrs 

5000 hrs 

Rotary Wing 

0 hrs 

0 hrs 

60 hrs 2500 hrs 0 hrs 

1000 hrs 

t 1 

Simulation Experience 







Fixed Based Simulators 

1 25 hrs 

2 5 hrs 

50 hrs- 125 hrs 

120 hrs 

250 hrs 

Moving Based Simulators 

I 

0 hrs 
J 

0 hrs 

1 

30 hrs 25 hrs 

130 hrs 

50 hrs 









The first day of running the subject generally reported 
after lunch. At that time, a standard description of the program 
objective, the mission and the landing modes was given the sub- 
ject. Next, the subject was seated in the cockpit and familiarized 
with the instruments with respect to the mission. Then, to be 
sure, the mission was understood a f ully-automatic mission was 
flown while another verbal description of the mission was made. 

At this point, the subject began a standard training sequence. 

The training runs began with the easiest mode to fly, and 
progressed to the most difficult. The mode order was: 

(1) Velocity Command with FDI 

(2) Velocity Command without FDI 

(3) Attitude Command II with FDI 

(4) Attitude Command II without FDI 

(5) Attitude Command I with FDI 

(6) Attitude Command I without FDI 

For each mode but the last, the subject was considered 
qualified on that mode when either two successful runs in a row 
were made or eight runs were made. The second provision of 
eight runs was made due to the time constraint. In general, 
pilots were only available for an afternoon and the next day. 

When the subject had qualified, two additional runs were made 
and then the subject went on to the next mode. For the last and 
most difficult mode only two out of four in a row were required 
for qualification. Subjects were given the option of a ten- 
minute break between each mode. The training sequence was generally 
completed within four hours. 

The production runs were made on the second day. For the 
production runs, the order of the modes was different for each 
subject to average out the effect of order on the data. Each 
mode was run until the subject requalified and then production 
runs were made. Five runs were made for each of the FDI modes 
and eight runs were made for modes without the FDI. The number 
of runs selected was based on the time constraint. The five/ 
eight ratio was chosen since the variance of the FDI runs was 
much smaller than those without FDI and therefore more runs without 
FDI were desired. Subjects were given the option of a ten-minute 
break between modes and lunch was taken between the third and 
fourth mode. The production sequence was generally completed 
within seven hours including an hour for lunch. 
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In addition to taking the performance data, the pilot was 
also asked the following question after finishing each mode 
(i.e., the mode just flown). "Would you fly this system for the 
defined mission in an operational environment?" The answers were 
to be simply, yes, no, or marginal. 


RESULTS 

The results of the experiment are shown for each subject in 
Figures 80 through 85. These data represent only those runs 
made on the second day after the training had been completed. 

In each figure, the overall performance index is shown for each 
system. Each individual data point is plotted and beside each the 
mean is circled and the standard deviation indicated. Below 
each plot is a table indicating for each system the number of runs 
made, the number of failures and, for the pilot subjects, the 
answer to the question: "Would you fly this system for the 

defined mission in an operational environment?" 

The composite results are shown in Figure 86. The means of 
each subjects performance index are plotted and beside them the 
composite mean is circled and the standard deviation of the 
individual means is indicated. Below the plot is a table in- 
dicating for each system the total number of runs made and the 
total number of failures. A summary of the pilot's answers to the 
question cited above is also included. 

As was expected, the relative performance of the automatic 
system is best. In the absence of electronic failure and sensing 
problems, an automatic system can be easily designed with only 
bandwidth considerations. As the level of pilot participation 
increases, the performance of the system decreases. This is due 
to the bandwidth limitations of the pilot. With the FDI, however, 
this decrease in performance is significantly reduced. This 
indicates that the FDI information source permits the pilot to 
operate at higher bandwidth than raw data sources. 

In addition to relative performance, an attempt was made to 
judge absolute acceptability of each system based on the percentage 
of approaches that failed to meet the minimum performance speci- 
fications. To be considered an acceptable operational system 90 
percent of the approaches made with a system had to be successful. 
The 90 percent figure was chosen based on current FAA Category-II 
system certification requirements as cited in reference 13. With 
this criterion the Attitude-I control system was unacceptable both 
with and without the FDI. All other systems were acceptable. These 
results correlated well with the pilots rating of the modes. All 
pilots rated both Attitude-I systems unacceptable or marginal while 
rating all other systems acceptable. 
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Figure 80.- Results of the experiment-Subject 1 
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Figure 82.- Results of the experiment-Subject 3 
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Figure 83.- Results of the experiment-Subject 4 
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Figure 84.- Results of the experiment-Subject 5 


105 



POINT AT 0.93 



J 


L VEL ATTI 
AUTO\ ATTH / 

v 

WITH FDI 


NUMBER OF RUNS 

8 

8 

8 

5 

5 

5 

5 

NUMBER OF FAILS 

8 

4 

0 


0 

0 

0 

ANSWER TO 
QUESTION 

NO 

YES 

YES 


YES 

YES 

NO 


Figure 85.- Results of the experiment-Subject 6 
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Figure 86.- Composite results of the experiment 
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In investigating the reason for the failure of both 
Attitude-I systems, the pilot workload must be considered. In 
the Attitude-I systems, the general comment of the pilots was that 
the workload was high. Since this was not the case with Attitude II 
and, since the Attitude-II system was acceptable, it can be de- 
duced that the prime cause for the high workload of Attitude I 
was the unaugmented vertical axis. The coupling of pitch and 
collective inputs, the nonlinear collective trim (i.e., power 
required curve) , and the effect of vertical gusts without aug- 
mentation made the vertical axis of Attitude I very difficult 
to fly. This is dramatically evident in the number of failures 
sustained in Attitude I with raw information sources. The effect 
of the FDI was to assist the pilot in flying the vertical axis 
but the workload remained high and the pilot was still unable 
to consistently remain within specified acceptable performance. 

In addition, the subject was forced to "tunnel" on the FDI (i.e., 
concentrate solely on the FDI) for this system and act solely as 
a servo. The high workload and the associated "tunneling" on the 
FDI led to the low pilot rating on this system. 

In considering the standard deviation of the various runs 
associated with each subject, it is seen that, in general, the 
lower the bandwidth required of the subject, the smaller the 
standard deviation. This is because as less action is required 
of the subject, less experience with each system is required to 
produce a consistent set of runs. In considering the standard 
deviation of the means presented in the composite, again, the 
lower the bandwidth required of the subject the smaller the 
standard deviation. This is because as less action is required 
of the subject, less general experience and innate talent is 
required to produce a consistent set of runs. These results on 
standard deviation are important with respect to training. They 
indicate that the lower the bandwidth required of the subject by 
a system, the lower the amount of training required to obtain 
proficiency with that system. 

Finally, each system was considered for the localizer track 
and ILS beam track phases of the mission alone. Such a mission 
would assume visual breakout while the helicopter is on the ILS 
beam at 150 ft of altitude. For this reduced mission, every run 
for every mode met the minimum performance specifications (i.e., 
all failures in the data occurred in the flare and deceleration 
maneuver). Therefore, all modes are considered acceptable 
operationally to 150 ft of altitude. This was to be expected, 
since such missions are currently flown with simple attitude rate 
damper systems . 
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CONCLUSIONS 


A fixed-base simulation evaluation of various low-visibility 
landing systems for helicopters has been made. The objective 
was to provide performance data for use in making trade-off 
studies required in selecting an operational system. Based upon 
the simulation, the following conclusions are made: 

(1) All systems considered are operationally acceptable 
for a straight-in 6-degree glide slope and localizer 
mission with visual breakout prior to 150 ft of altitude. 

(2) For the total low-visibility mission with visual break- 
out prior to 100 ft of range, the attitude command mode 
with an unaugmented vertical axis (i.e.. Attitude Mode 
I) is operationally unacceptable with raw information 
sources or with a flight director indicator (FDI ) . All 
other systems are operationally acceptable for the total 
low-visibility mission. 

(3) For the total low-visibility mission, augmentation of 
the vertical axis to provide a tight Rate of Descent 
Command system is essential for an acceptable system. 

(4) There is a substantial improvement in performance of the 
Velocity Command Mode over the attitude mode with 
vertical axis augmentation (i.e.. Attitude Mode II). 

(5) The Automatic System has the best performance. 

(6) An FDI improves the performance of each of the three 
control modes. This is done, however, at the expense 
of "tunneling" on the FDI with resultant loss of 
knowledge of system status. 

(7) The lower the bandwidth required of the subject by a 
mode, the less the amount of training required to obtain 
proficiency with that mode. 
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APPENDIX A 


The flight control modes were defined by considering the 
complexity of the sensors and computations which would be required 
to implement each mode in operation. The sensors and associated 
filters were only hypothesized. They were not studied in depth 
nor were they simulated. 

First, three sets of basic sensors and associated instru- 
ments were defined as follows: 

VFR Flight Sensors and Instruments 
Barometric Altimeter (Displayed) 

Barometric Rate of Descent (Displayed) 

Indicated Airspeed (Displayed) 

Three Angular Rate Gyros (for Stabilization) 

Sideslip Sensor (for Stabilization) 


IFR Flight Sensors and Instruments 
Vertical Gyro (Displayed via Attitude Indicator) 
Directional Gyro (Displayed) 


Landing Sensors and Instruments 
Glide Path and Localizer (Displayed) 

Range and Range Rate (Displayed) 

Radar Altimeter (Displayed) 

From these sensors the flight control modes were defined 
as follows. 

Attitude Mode I - Would use the information from the VFR 
and the IFR sensors. 

Attitude Mode II - Would use the information from the VFR 
sensors, the IFR sensors, and, in addition, require at 
least a set of three body-mounted accelerometers. The 
accelerometer data would be used to blend with the 
barometric altimeter to derive a good vertical velocity 
estimate. A hypothesized filter is shown in Figure A-l. 
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Figure A-l.- Hypothesized vertical velocity filter 


Velocity Mode - Would use the information from the VFR 
sensors, the IFR sensors, the accelerometers, and the 
landing sensors. The accelerometer data would be used 
to blend with the landing sources to derive good V x and 
Vy information. The filters could be similar to that 
for the vertical velocity. 

Automatic Mode - Would use the same information as the 
velocity mode. 
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APPENDIX B 


The stability and control derivatives for Flight Condition 1 
are listed in Tables B-I and B-II for the longitudinal and 
lateral-directional modes, respectively. The operating conditions 
are : 


Gross Weight 

13,400 lbs 


Center-of-Gravity Position 

Normal 


I 

9,203 slug 

ft 2 

XX 


I 

75,914 slug 

ft 2 

yy 


i 

zz 

71,786 slug 

ft 2 

J 

-7,114 slug 

Hi 

ft 

K) 

xz 


Rate of Descent 

0 ft/mi n 


Altitude 

Sea Level 
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TABLE B— I . — LONGITUDINAL STABILITY DERIVATIVES 
AT FLIGHT CONDITION I 


Forward 

Velocity 

kts 

0 

20 

40 

60 

80 

100 

120 

140 

Forward 

Velocity 

ft/sec 

0.0 

33.76 

67.52 

101.28 

135.04 

168.80 

202.56 

236.32 

a TR 

deg 

0.0 

8.13709 

6.54894 

4.67298 

2.31058 

1.40310 

-2.27505 

-6.6963 

0 TR 

deg 

9.30627 

8.19834 

6.62235 

4.75227 

2.32294 

1.43387 

-2.26671 

-6.6331 

V m 

ft/ sec 2 
ft/sec 

-.02540 

-.00181 

-.02156 

-.03604 

-.04642 

-.05579 

-.06456 

-.07206 

V ra 

ft/sec 2 
ft/ sec 

.05449 

.06818 

.08255 

.08944 

.08548 

.10343 

.08523 

.05666 

X q /m 

f t/sec^ 
rad/sec 

.60185 

.74915 

.87508 

.84957 

.73079 

1.24872 

.77601 

.25183 

X 6e /m 

ft/sec 2 

in 

.17696 

.13988 

.12312 

.14237 

.16406 

-.05165 

-.03776 

-.01954 

W" 

ft/sec 2 

in 

1.20482 

.97467 

. 87948 

.80253 

.68335 

.88200 

.72040 

.48341 

Z u /m 

ft/sec 2 

ft/sec 

.06009 

-.12594 

-.08296 

-.02192 

.01396 

.06077 

.06546 

.04980 

V ra 

ft/sec 2 

ft/sec 

-.36933 

-.48399 

-.63639 

-.80152 

-.92055 

-1.00063 

-1.04552 

-1.10516 

V m 

f t/sec^ 
rad/sec 

-.71511 

-1.16872 

-1.77844 

-1.81400 

-1.81986 

-2.22039 

-2.25955 

-2.47935 

Z 6e /m 

f t/sec^ 
in 

-.00407 

.21188 

. 51943 

. 56820 

.52527 

.46341 

.41480 

.36743 

Z Sc /m 

_ / 2 
ft/sec 

in 

-7.43006 

-7.23138 

-7.65410 

-8.52446 

-9.49005 

-10.2660 

-11.0874 

-11.7111 

M u /I yy 

2 

rad/ sec 
ft/sec 

.00656 

.00645 

-.00587 

-.00670 

-.00582 

-.00185 

- . 00120 

-.00089 

M v/ I yy 

2 

rad/sec 

ft/sec 

-.00285 

.00978 

.01630 

.01363 

.01154 

.00956 

.00774 

.00691 

V I yy 

2 

rad/sec 

rad/sec 

-.73173 

-.96002 

-1.31158 

-1.45996 

-1.52219 

-1.62003 

-1.59067 

-1.51570 

M 6e /I yy 

2 

rad/ sec 
in 

. 35447 

.35364 

.40144 

.45022 

.48135 

.51572 

.53335 

.54565 

M 5c^ I yy 

2 

rad/sec 

in 

-.04765 

-.04252 

.04556 

.06776 

.06725 

.04707 

.03948 

.03505 

e TR 

in 

.66523 

-.06503 

-.23516 

.28888 

.77817 

-.06119 

.16392 

.34298 

6c tr 

in 

5.01959 

4.47346 

3.73135 

3.51111 

3.84917 

4.67777 

6.04983 

8.02025 
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TABLE B-II 


LATERAL STABILITY DERIVATIVES AT FLIGHT CONDITION I 


Forward 

Velocity 

kts 

0 

20 

40 

60 

80 

100 

120 

140 

Forward 

Velocity 

ft/sec 

0.0 

33.76 

67.52 

101.28 

135.04 

168.80 

202.56 

236.32 

V m 

ft/sec 2 

ft/sec 

-.026638 

-.05408 

-.07986 

-.11823 

-.15121 

-.18040 

-.19766 

-.23442 

Yp/m 

ft/s ec^ 
rad/sec 

-.76514 

-.95618 

-1.15981 

-1.24824 

-1.23341 

-1.02295 

-.76619 

-.29372 

Y r /m 

ft/sec 2 
rad/ sec 

-.12517 

-.18179 

-.20082 

-.12813 

-.05819 

.13380 

.12514 

.22125 

Y «a /ra 

ft/sec^ s 
in 

.99794 

.99700 

.97673 

.96426 

.96955 

.97123 

1.01613 

1.10225 

Y 6r /ra 

ft/sec 2 

in 

.14652 

.13634 

.11652 

.10583 

.09849 

.11641 

.13883 

.15262 

V^x 

2 

rad/sec 

ft/sec 

-.00778 

-.01305 

-.01388 

-.01720 

-.02442 

-.03188 

-.04724 

-.05080 

L p /I xx 

rad/sec 2 

rad/sec 

-.50730 

-.57483 

-.64045 

-.65271 

-.62549 

-.52243 

-.41134 

-.23476 

V 1 ** 

2 

rad/sec 

rad/sec 

-.02297 

-.04571 

-.05943 

-.02270 

.01172 

.09531 

.11659 

.19955 

L 6a /:E xx 

rad/sec 2 

.46536 

.46558 

.45954 

.45575 

.45805 

.45820 

.47354 

.50528 

in 

L «r/ J xx 

rad/sec 2 

in 

-.12638 

-.13036 

-.13422 

-.13544 

-.13910 

-.13256 

-.13248 

-.14433 

V Z z 2 

2 

rad/sec 

ft/sec 

.00013 

.00037 

-.00123 

-.00290 

-.00415 

-.00536 

-.00377 

-.00060 


2 

rad/ sec 
rad/ sec 

-.01831 

-.02076 

-.02598 

-.03947 

-.05344 

-.06864 

-.07797 

-.07471 

N /I 
r zz 

2 

rad/sec 

-.05847 

-.05450 

-.04597 

-.05020 

-.05438 

-.05078 

-.08642 

-.12617 

rad/ sec 

"i.^zz 

2 

rad/sec 

in 

.03001 

.02923 

.02764 

.02663 

.02633 

.02696 

.02857 

.02970 

« , «r /1 z* 

2 

rad/ sec 
in 

.17584 

.17578 

.17261 

.17032 

.17141 

.17146 

.17933 

.19536 

6 

a TR 

in 

.12983 

.08462 

.09191 

.13466 

.17931 

.29754 

.44154 

.49785 

6 r 

TR 

in 

-.17764 

-.04701 

-.08508 

-.33705 

-.61293 

-.75397 

-1.04713 

-1.14123 
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APPENDIX C 


The formal presentation of the wind model follows. Given 
the differential equation: 


f = -3 r + cr $ /2 ? n 

n n n 


(Cl) 


where 3 n > 0 and £ n is white noise with a spectral density equal 
to one. It can easily be shown that the steady-state auto- 
correlation function is: 


E(r(t)r(t+x)) 


2 

a 

n 



(C2 ) 


If one wishes T to have a variance, P , and a correlation time, 
a n , then: " ! 


6 n " Va n 


a n " /a n p n 


(C 3 ) 


To approximate this process on a digital computer, Eq . (Cl) was 
time discretized and a white, gaussian-distributed sequence 
[C™] was generated to provide a piecewise constant approximation 
to £ n (t). Under the assumption that the product of the dis- 
cretization interval, A n , and of 3 n is small compared to unity, 
Eq. (Cl) becomes: 


r n+1 = (1-A 3 ) T n + a g /2 A C n 
n n n n n n 


(C4 ) 


where superscripts indicate time instants. Since C n (t) is the 
formal derivative of a Wiener process, its integral must possess 
the same statistical properties as the Wiener process, that is: 


e 


e 


(C c " <o) d °) ■ ° 

(Ef c " (a) da ] 2 ) = lk 


(C5 ) 
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The sequence 


u n ] 


must be chosen such that 


e 



rt/A t/A 

EE" 

Li=l j-1 



t 


(C6) 


Since [c“] is white, equation (C6) becomes 



var (? 2 ) A n | t | 


(C7) 


Upon comparing Eq. (C7) with Eq . 


(C6) , it is found that 


VAR (sjj) = 1/A n 


(C8) 


is required to provide the desired noise sequence. Defining 
[£ n ] = [n n //A n ] where [n£] is a white, gaussianly distributed 
random sequence with zero mean and unit variance, causes Eq. (C4) 
to become 


1-11+1 - U-Vn’ 1- " + °n B n 


/2A~ r] n 
n n 


(C9) 


The sequence (ri£) is generated by summing twelve samples from a 
normalized, uniform distribution in (0, 1) and then by subtracting 
off the mean value of the sum. The Central Limit Theorem states 
that the distribution of these sums will be approximately 
gaussian. The uniform distribution was generated by using a 
multiplicative, congruential random number generator tailored 
to the SDS 9300, which has a word length of 24 bits. The 
[n£] is given by 

5 1 ? , = 2899* H n (mod 2 23 ) 
l+l i 

... — 0 _ 

with = 7 

n n = if 2 -23 “5 1 - 6 n = 0,1,2,--- 

i=l 
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